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SUMMARY 

This report documents the aerodynamic and mechanical design of 
a radial inflow turbine stage for a Brayton rotating unit application. 

The work was performed under NASA Contract No. NAS3-13492. The new tur- 
bine is an alternative to one previously designed and tested. 

Since the single-stage turbine in a turbo-alternator-compressor 
unit supplies both alternator and compressor drive power and the power re- 
quirements of the latter are significantly greater than that of the former, 
a one-percentage-point change in turbine efficiency will produce a rela- 
tively large change in the useful power output of the unit. In this par- 
ticular turbine application, the combination of stress levels and operat- 
ing temperature are such that the primary design emphasis could be placed 
on achieving an aerodynamic optimum within the over-all constraint of the 
life requirement of the turbine. 

The work reported consists of three phases; the preliminary 
aerodynamic and mechanical design, the final aerodynamic design, and the 
final mechanical design including the production of manufacturing draw- 
i ngs. 

The selected stage configuration differs from the previously 
tested stage in a number of respects. The rotor design features back- 
slope blading at rotor inlet and a tandem blade row design, the meridional 
flow path through the nozzle row is flared, the downstream diffuser sec- 
tion was designed for a linear static pressure rise, and the inlet scroll 
over-all dimensions were reduced compared with the original. Compared 
with the existing design, the efficiency of the stage is predicted to be 
between 2 and 4 percentage points higher at the design-point over-all 
tota 1 -to-stat i c pressure ratio. The relatively large band of uncertainty 
assigned to the predicted efficiency is directly associated with the cur- 
rent deficiencies of loss prediction procedures for radial inflow tur- 
bines in general and their application to the nozzle row and rotor con- 
figurations of the new design. 

The mechanical design of the rotor was undertaken for the BRU 
application. However, the stationary parts were designed for cold flow 
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testing on the NASA test stand. Drawings for all new parts for the tur- 
bine stage test have been reduced and reproduced for inclusion in this 
design report. 
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INTRODUCTION 


Problem Statement 

The basic objective of this program was to design a radial in- 
flow turbine that surpasses the performance of an existing unit (Refs 
1 and 2) and which could replace the original with a minimum of modifica- 
tions both in an actual BRU and the NASA turbine test stand. 

The turbine drives a centrifugal compressor and an alternator 
in a closed Xe-He Brayton cycle. Since the compressor absorbs most of 
the turbine power output., a one-point improvement in the turbine effi- 
ciency will produce an approximate 3 per cent improvement in cycle ef- 
ficiency. The principal objective was to achieve the highest possible 
stage efficiency for the over-all design-point tota 1 -to-stati c pressure 
ratio; the stage comprises an inlet scroll, a nozzle row, the rotor, and 
an exhaust diffuser. 


Design Requirements 

The design requirements for the turbine are as follows: 


Worki ng F 1 ul d , Hel i um- 
Xenon Mixture 

Turbi ne I nlet Total 
T emperature 

Turbine Inlet Total Pressure 

Specific Heat at Constant 
Pressure 

Ratio of Specific Heats 
Molecular Wei ght 
Turbine Rotative Speed 
Direction of Rotation 

Wei ght Flow 

Turbine Inlet Total-to- 
Di f fuser Exi t Stat i c 
Pressure Ratio 


2060 deg R (1144 deg K) 

25.0 psia (172 kN/m 2 abs) 

0.05925, Btu/lb/deg R 
(248 J/kg deg K) 

1 . 667 
83.8 

36.000 rpm 

Clockwise Looking Up- 
stream at Rotor Exit 

0,7484 Ibm/sec 
(0.3395 kg/sec) 

1 .763 



Turbine Inlet Total-to- 
Diffuser Exit Total 
Pressure Ratio 


1.7^9 or Greater 


Technical Approach 

The over-all technical approach to this problem was outlined by 
NASA personnel initially and further developed by NREC during the prellmi 
nary design effort. Considerable flexibility in the rotor design was per 
mitted by the moderate wheel speed and the relatively low operating tem- 
perature of the closed-cycle rare gas turbine. Suggested rotor design 
modifications impl i ed performance gains in both the vaned nozzle and the 
di f fuser. 

For the rotor the initially recommended design change was an 
increase in outer diameter over the previous design, which at constant 
rotative speed corresponds to increased tip speed. The consequent re- 
duction in the absolute tangential inlet velocity would result in the use 
of backward swept blades which, in turn, would permit a more uniform 
loading distribution along the blade surfaces when compared to radial 
blades. A reduction in the number of rotor blades was also considered 
desirable in order to obtain the optimum balance between friction and 
diffusion losses. For this reason, alternate blading concepts to the 
conventional splitter configuration were also considered worthy of evalua 
tion. In addition, the other dimensions of the rotor such as inlet pass- 
age width, passage length, and outer diameters at inlet and exit were 
considered variables in the design optimization. 

The nozzle ring would benefit from an increase in discharge 
radius corresponding to that of the rotor, due to the reduced discharge 
velocity level. Here again, the remaining major nozzle dimensions (in- 
let radius, passage width) and vane number were to be selected in such a 
manner that the combined effect of loading, friction, and secondary flow 
losses would be minimized. Initial NASA design recommendations were 
shorter vane chords and flared nozzle walls. 

Significant performance improvement was expected from a rede- 
signed diffuser, since the experimental evaluation of the previous com- 
ponent (Ref 1) had indicated substandard performance. Diffuser 
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effectiveness as a function of diffusion is clearly related to mean rotor 
discharge conditions and as such are involved in the stage design opti- 
mization. Beyond that, however, the detailed design of the diffuser 
should include concepts that would reduce the adverse effects of the un- 
settled flow field emerging from the rotor. 

Report Arrangement 

The systematic design of the subject turbine evolved from pre- 
liminary aerodynamic and mechanical analyses to the detailed aerodynamic 
design and mechanical analyses. The report follows this order in the 
main part. There are three appendices; the first and second document 
the design and analysis procedures for the aerodynamic and mechanical 
design, respectively, while the third contains reproductions of the draw- 
ing reduced in size for report purposes. 
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PRELIMINARY DESIGN ANALYSIS 


Objectives and Approach 

The purpose of the preliminary analysis is to determine the 
principal dimensions and flow conditions (vector diagrams) at the inlet 
and discharge stations of each component as shown in Figure 1. The 
evaluation must consider aerodynamic performance and mechanical feas- 
ibility as a function of these dimensions and flow conditions and es- 
tablish a basis for their final selection, the relevant criterion being 
optimum stage efficiency within specified mechanical constraints. 

To achieve this purpose it Is necessary to identify the sig- 
nificant geometric variables, to define the relevant stress criteria, 
and to formulate the relationship between geometry, flow field, and 
pressure losses on the one hand, and that between geometry and stresses 
on the other. The aerodynamic relationships were established by means 
of simplified flow field representation and suitable expressions for 
distinct loss mechanisms associated with friction, diffusion, dissipa- 
tion, and nonuniformity of the flow field (loading), henceforth summarily 
termed the loss system. The mechanical relationships were similarly 
found from simplified expressions for rotor disk and blade root bending 
stresses and for blade natural frequency and vibratory stresses. Hence, 
a parametric analysis of performance and stress could be used to provide 
a rational approach to a preliminary design selection. 

Preliminary Aerodynamic Analysis Procedure 

A review of recent literature related to radial inflow tur- 
bines showed that the main contributions to radial turbine technology 
since the completion of NREC’s major turbine program (Refs 3 and 4) are 
the reports on the turbine which was being redesigned and the various 
reports on internal blade row flow analysis procedures of Katsanis (for 
example, Ref 5)* No improved loss correlations have been recently de- 
veloped, and data on which to base improved correlations are very limited. 
Further analysis was therefore based on the methods presented in Refer- 
ences 3 and 4. 
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Flow Analysis Procedure 

The flow analysis consists principally of satisfying continuity 
and work requirements for the given design flow, pressure ratio, and 
specified geometry, taking into account the losses that occur across 
each component. Since the losses depend on the flow quantities to be 
calculated, the solution of the final flow quantities is arrived at by 
an iterative procedure. The iteration terminates when the change in ef- 
ficiency from one iteration to the next iteration is within a prescribed 
limit. For the first iteration quantities which are not known are either 
assumed or approximated. In the succeeding iterations, the unknown quan- 
tities are taken to be those calculated In the previous iteration. 

As part of the stage efficiency calculation, continuity is sat- 
isfied at inlet and exit stations of the rotor and nozzle vanes. At these 
stations either the tangential velocity or flow angle is known. From the 
given area, design flow and the estimated total pressure and temperature 
continuity is solved iteratively to calculate the flow quantities. 

As a first step in the efficiency iteration loop, the exit total 
temperature is calculated from the specified pressure ratio and the ef- 
ficiency calculated In the previous iteration. The diffuser exit static 
pressure is obtained from the given inlet tota 1 -to-exi t static pressure 
ratio and inlet total pressure. Since the design flow is specified and 
the annulus area is known from the dimensions, the diffuser exit flow 
conditions are calculated assuming no discharge swirl. 

The diffuser inlet total pressure Is calculated from the diffu- 
ser loss formulation for the known area ratio, prescribed velocity gradi- 
ent at the diffuser inlet, and che inlet blockage estimated from the aver- 
age wall momentum thickness at the rotor exit. Since the absolute flow 
angle at the rotor exit is specified, the flow conditions are obtained by 
continuity iterations. 

The rotor losses are calculated as the sum of the clearance 
loss and the combined friction and diffusion losses. The rotor inlet 
total pressure is obtained from the work and the calculated losses. The 
absolute inlet tangential velocity is calculated from the work, the exit 



8 


tangential velocity, and the blade speed, The remaining rotor inlet flow 
conditions are then calculated by continuity iterations. 

The nozzle vane exit total pressure and tangential velocity are 
calculated from the vaneless space loss and conservation of momentum. The 
continuity iteration is then employed to calculate the remaining flow con- 
ditions at the exit of the nozzle vanes. The nozzle vane inlet total pres 
sure is calculated from nozzle vane losses. Since the absolute flow angle 
at the nozzle vane inlet is prescribed, the flow conditions are obtained 
by continuity iterations. 

The scroll or the turbine inlet total pressure is then obtained 
from the calculated scroll losses. The turbine efficiency is estimated 
based on the calculated inlet total pressure. The efficiency is then com- 
pared with the efficiency of the previous iteration. If the efficiency 
difference is greater than certain prescribed amounts, the new efficiency 
is assumed and the iteration procedure repeated with the given pressure 
ratio until the efficiency difference falls within the prescribed limito 

Loss System 

During an initial evaluation of alternative loss systems, it was 
shown that the axial turbine loss correlation of Reference 6 could be modi 
fied to accurately predict the design-point efficiency of the existing BRU 
turbine as reported in Reference 1. However, while the over-all perfor- 
mance was accurately predicted, the individual component losses were less 
accurately predicted. Hence, although the initially proposed correlation 
is well suited to a parametric analysis, it was considered necessary to 
further refine the loss system in order that valid predictions of modifi- 
cation to the existing design geometry would be obtained. The measured 
static pressure in the stator-rotor gap was used to obtain estimates of 
the stator and rotor performance. Even though the calculation of a stator 
loss coefficient by means of static pressure data is subject to possible 
inaccuracy, it was clear that neither the axi a 1 -turb i ne -based loss corre- 
lation of Reference 6 (NREC Report 1125-1) nor the radial turbine stator 
correlation of Reference 3 (NREC Report 1067-1) would predict the stator 
loss of the existing turbine; the "1067" prediction gave too low a loss, 
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and the "1125" prediction gave too high a value. A more recent correla- 
tion, obtained from planar cascade tests of low-aspect-ratio axial turbine 
blade sections, was therefore considered. This correlation is based on 
an NREC program for USAAVLABS reported in Reference 7 (NREC Report 1137 _ 0 
(The report will be issued as AVLABS Report No. 70— 1 . ) In applying this 
third correlation, the stator system loss was considered to have three com 
ponents, a scroll loss, the stator blade row loss, and the vaneless space 
loss. From the following tabulation it will be seen that the 1 1 1 137" and 
"1067" prediction produces the best agreement with the value derived from 
test data. 


Loss Coefficient Derived 
from Experimental Data 


Predicted Loss Coefficients 


= 0.082 
"1067" = 0.036 
"1125" = 0. 148 
"1137" + "1067" 


0 . 122 


It appeared that a logical modification of the axial turbine 
blade row correlation to include the effect of the radius change through 
the stator blade row produced good agreement between the predicted total 
pressure loss and that derived from the experimental data of Reference 1. 

The selected nozzle row total pressure loss coefficient is as 

fol 1 ows : 


AP a _ C |ba.~ftO<i-si ^Q. 007 4- 0.044^j7 S 4 ^ \ 

P oa-P 2 ” 1 0.6+o.8C0S<*J l \£JL 

So. U4 + | < ' ) 


Equation 1 is a modified version of that obtained for axial 
turbine cascades in Reference 6. The modification introduced for the 


radial nozzle is the diameter ra 


. . ©2. . . 

fo ■=- ; it 


©1 


was reasoned that the con- 


servation of angular momentum in a radial space should be taken into ac- 
count when calculating the row loss due to tangential loading. The four 
remaining terms of Equation 1 introduce loss level correction for row 
acceleration, Vi k , passage aspect ratio, AR , trailing edge blockage, 
£e/s , and Reynolds number. 
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Introduction of the above loss formulation Into the radial tur- 
bine prediction program produced acceptably close agreement between test 
and prediction. The total pressure loss coefficient predicted for the 
scroll, nozzle vanes, and vaneless space is 0.0745 compared with the 0.082 
value derived from the experimental data. In view of the possible inac- 
curacy associated with deriving a total pressure loss coefficient from 
measured static, the agreement was considered satisfactory. 

In the case of the rotor losses, the initially selected loss 
correlations, based on Reference 3, considerably underestimated the rotor 
losses as deduced from the stage tests. At this point, a number of pos- 
sible explanations were sought rather than some arbitrarily selected "cor- 
rection 11 factor introduced. Two basic deficiencies of the original loss 
system were identified. One, the loss system did not take splitters into 
account other than by replacing blade and splitter numbers by an equiva- 
lent number of full blades. Two, since the loss formulation assumed a 
linear variation of relative velocity from inlet to exit, it could not 
account for a major loss associated with a deceleration along the hub 
line (which is inherent to most radial inflow turbines). Accordingly, a 
loss formulation was developed which introduced an intermediate station 
(Station 31, see Fig 1) between rotor Inlet and rotor exit (Stations 3 
and 4, respectively). Velocities calculated at the intermediate station 
are used in designs both with and without splitters. For the existing 
design, which employs splitters, two sets of velocity calculations are 
performed at Station 31 with and without the blockage of the splitter. 

The calculation station in which the splitter blockage is included Is 
named 331 (and is used to calculate losses from 3 to 31) • The station 
immediately downstream of the splitter is identified as Station 314 (and 
is used to calculate losses in the section between 31 and 4). 

The calculation is based on a combined friction and diffusion 
loss. Velocity distribution calculations are used to calculate diffusion 
ratios which are then related to boundary-layer displacement thickness 
and, hence, a rotor total pressure loss coefficient. For the purpose of 
the calculation, the rotor flow is divided into inner and outer flows. 
Conditions for the inner flow are assumed to be the average of the mean 



line and hub flow conditions. Similarly, the outer portion of the flow 
is represented by the average of the mean and outer streamline flow con- 
ditions. For the full blade, four components of diffusion are identified, 
that is, 


n <1- §£!*£)+ (± - W 3i^ )+fl— V/i- 

w v w 3,« y ^ w S3lj y v. w i 


w. 




where W 3’ W 33l’ W 3I4' and W 4 are relative velocities at the four cal- 
culation stations. The subscripts U) and denote wall and average val- 
ues, respectively. Four values of correspond! ng to the inner and 
outer, pressure, and suction surfaces are calculated. Only relative ve- 
locity ratios less than unity are used in the calculation. Each of the 
four values of D w is used to calculate a displacement thickness, , 
where 


°z — — I .a < D ^ 0.4- 

S = 

1 (1 -D^- 5 


D > O. 4 

w 


(3) 


The four (Rvalues are used to obtain an average value , which in 

turn is used to calculate the loss coefficient, ,AP /D _ P > given 

L _ ° f r O*c r 4 

L 


by 


AR. 


= 4 K 


Dr 


w 


P04 P(. 

where i s a secondary flow factor and is the hydraulic diameter for 

the rotor exit flow area. 

The above loss formulation was used in the parametric analysis. 


Preliminary Mechanical Analysis Procedure 


Procedures used in the mechanical phase of the program are 
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fully documented in Appendix II. The following topics represent areas of 
primary concern in the preliminary design analysis of alternative rotors: 

1. Blade root bending and tensile stresses. 

2. Centr i f uga 1 1 y loaded disk stresses. 

3. Blade natural frequencies. 

Each of these topics is briefly discussed below. 

Blade Root Stresses 

Experience has shown that where blades with any significant degree 
of backslope are used, the blade root bending moment at the disk periphery 
represents the stress limiting situation. If there is lean on the axial 
portion of the blade, however, the bending stresses produced in this sec- 
tion could very quickly become limiting. Tensile stresses are very rarely 
1 imi ti ng. 

Values of stress, when obtained, are compared with the appro- 
priate strength value, usually yield strength, of the material to be used. 

Disk Centrifugal Stresses 

In many instances, the disk Is not the limiting item in stress 
considerations. In addition, considerable improvement in both stress and 
deformation patterns may be affected by proper shaping of the impeller 
backface. The equations used In the preliminary analysis consider the 
two important stresses. These are the basic disk maximum stress and the 
stress induced by the moment exerted on the disk by backward sloped 
blades. The former was generally the more important, and it is deter- 
mined entirely by the impeller tip speed and the material properties. 

A conservative solution is represented by a solid homogeneous 
circular disk of uniform thickness. For a solid circular disk or cylinder, 
the maximum radial and tangential stress occurs at the center of the disk. 
Formulas for stress calculations are those of Reference 8. 

Blade Natural Frequencies 

Blade natural frequencies were predicted using the methods and 
equations of Reference 9. The primary consideration in the preliminary 
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design stage was that any blade design likely to be selected as an aero- 
dynamic optimum would not have a calculated natural frequency less than 
1.3 times the highest fundamental frequency. The highest fundamental 
frequency was assumed to be the number of rotor or stator blades, de- 
pendent on which was the highest number, times the design rotational speed 


Parametric Investigation 


The parametric study consisted of selecting three datum designs 
and systematically varying the previously identified design geometric 
parameters. The latter datum designs were selected from the parametric 
investigation of the previous datum designs and only the parameters which 
showed improvement in turbine efficiency were used for parametric analy- 
sis of the latter datum designs. The major parameters of the three datum 
designs are listed in Table I. 

Eleven design parameters were initially identified for parametric 
investigation. They consist of five diameters, three passage widths, a 
number of nozzle vanes, a number of rotor blades, and the axial length of 
the rotor. The five diameters correspond to the nozzle vane inlet and 
exit, rotor inlet, and rotor exit hub and shroud. The three widths are 
at the nozzle vane inlet and exit and rotor inlet. For the main portion 
of the parametric study, the number of independent variables was reduced 
to nine. 


The vaneless space is provided mainly for the vane wakes to 
mix out before entering the rotor. Any large space would introduce ad- 
ditional frictional loss. Since the diameter ratio (vane exit to rotor 
inlet) for wel 1 -des i gned nozzle vanes is approximately 1 . 04 (Ref 3) and 
the reference turbine (Ref 1) has a diameter ratio of 1.03, the diameter 
ratio for the parametric study was assumed to be 1.03 for the three datum 

designs". Hence, the nozzle exit diameter, £l , was eliminated as an 

2 . * 

independent variable for the initial parametric analysis . Since the 


ic 


The influence of the ratio Vpj on stage performance was 
for the final selection of stage geometric variables, 


I nvesti gated 
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length of the vaneless space is small compared to the length of the noz- 
zle vanes, the nozzle vane exit passage width, b^ , was assumed to be the 
same as rotor Inlet width, b^. Hence, the width b^was removed from 
the independent variables. 

The first datum turbine (Datum 1) was selected with a 3 per cent 
increase in nozzle vane inlet diameter, , and a 4 per cent increase in 
both rotor inlet diameter, X) g > and nozzle vane inlet width, , com- 
pared with the reference turbine- Datum 1 turbine was investigated for 
designs both with and without splitters. For the case without the split- 
ters, the number of blades was taken to be the equivalent number based 
on meridional length. The efficiency of the Datum 1 turbine without 
the splitter was found to be about 0.6 per cent less than that of Datum 1 
with the splitters. Therefore, splitters were assumed in all designs 
during the major portion of the parametric analysis, 

in the parametric study all the variables except the rotor exit 
hub diameter and the number of blades were varied +4 per cent from the 
values of the datum designs. The rotor exit hub diameter was increased 
with the shroud diameter such that the exit area remained constant (the 
product of the mean diameter and the blade height being held constant). 

The results of the parametric investigation are shown in Fig- 
ures 2 through 6. In all but one case the variation in over-all total- 
to-static efficiency is shown as a function of the particular parameter. 
Each figure presents results for the first datum design, while some also 
contain data for Datum Designs 2 or 2 and 3, depending upon the relative 
importance of the parameters. Since the stress data are principally de- 
pendent on the tip diameter (or rotor tip speed), these have only been 
included with the fluid dynamic efficiency results of Figure 2 where the 
results of varying]} as an independent parameter are shown. 

Considering the five figures in more detail, the effects of 
change in nozzle vane inlet diameter, , and rotor inlet d i ameter , 

are shown in Figure 2. An increase in Di increases the efficiency as a 
result of changes in reaction and the Reynolds number of the nozzle vanes. 
Similar trends are exhibited for Datum 1, 2, and 3 turbine designs, but 
the increase in of 0.2 per cent in the case of Datum 1 decreases to 
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0.16 per cent for Datum 3 when Pi is increased by 4 per cent from the 
datum values. This reflects the reduced level of stator loss for Datum 3 
turbine compared to the Datum 1 stator loss. The effect of change in 
is also shown in Figure 2. The efficiency drops when is decreased 

as a result of the predominant increase in rotor diffusion.. For an in- 
crease in , the decrease in rotor diffusion loss is small compared to 
the increase in rotor friction and nozzle vane loss, the latter of which 
is the result of higher nozzle vane exit diameter, Therefore, an 

increase or decrease of 4 per cent in reduces the turbine efficiency. 
The efficiency attains a maximum at a diameter ratio slightly under 1,02. 
The stress levels are within preset limits for this diameter ratio. There- 


fore, for the other two datum turbines the diameter, D 


3 s 


be 1.02 times the Datum 1 value and held constant. 


was chosen to 


The variation of efficiency with rotor exit shroud diameter, 

, as well as with the mean diameter, , for constant exit annulus 

area is shown in Figure J>. A decrease of 4 per cent, in decreases dif- 

fusion losses in rotor by approximately 15 per cent. The higher diffusion 
imposed on the diffuser results in a diffuser loss which is more than dou- 
ble that of the Datum 1 turbine. The net effect is a drop of 1.4 per cent 
in turbine efficiency. If, on the other hand, the diameter, De* , is in- 
creased, the increase in rotor diffusion loss predominates over the de- 
crease in diffuser loss and brings down the turbine efficiency by 0.3 per 
cent. The turbine efficiency is predicted to steadily drop when the rotor 
exit mean diameter is increased as a result of increased rotor diffusion 


loss. This effect is opposed by the beneficial effect due to reduced up- 
stream diffusion anticipated at the rotor hub. A review of the data shown 
in Figure 3 led to the conclusion that the trade-offs involved in changes 
of the exit dimensions of the rotor would only be established as a result 
of a detailed design analysis rather than by parametric analyses involving 
the geometries of the rotor exit. 

The change in turbine efficiency with variation of nozzle vane 
inlet width, , is presented in Figure 4. For all three datum turbines 
the efficiency increases with increase in reflecting the increase in 

reaction across the nozzle vanes. The effect of change in rotor inlet 
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passage width on ""^gg' 5 a ^ so included in Figure 4 . It can be seen from 

the figure that the efficiency increases with increase in , but the 

3 

relative change is small especially for the Datum 2 design. This increase 
is the net result of increased losses in nozzle vanes caused by lower re- 
action and decreased rotor losses as a result of higher reaction and lower 
average relative velocity levels. 

Figure 5 presents the effect of change in the number of both 
nozzle vanes and rotor blades. When the number of nozzle vanes is in- 
creased, the aspect ratio losses come down while the trailing edge losses 
go up, resulting in a net increase in \s . The rate of increase in ef- 
ficiency becomes negligible at higher number of nozzle vanes. in the case 
of rotor blades the efficiency decreases for both increase and decrease 
in the number of blades. When the number of rotor blades is increased or 
decreased, the decrease or increase in clearance loss approximately com- 
pensates the increase or decrease in friction loss. The total diffusion 
loss increases for lower numbers of blades; the total diffusion loss shows 
a slight increase even though the diffusion loss per blade is decreased. 

The efficiency decreases with increase in axial length, as shown 
in Figure 6. This decrease results from higher friction and diffusion 
losses. The clearance loss remains practically constant with rotor axial 
distance. The effect of vaneless space ratio (nozzle vane exit and rotor 
inlet) was investigated for the final selected design only. When the 
vaneless space diameter ratio ( was increased, the nozzle vane 

diameter ratio ( was k e Pt constant. As shown in Figure 6, the 

efficiency decreases with increase in as a result of increased 

losses in the vaneless space. The nozzle vane and scroll losses remain 
essentially constant. 


Preliminary Design Selection 

The following summary presents the selection of the principal 
design variables and the basis of the selection. With the exception of 
the rotor blade number the parameters selected following the parametric 
analysis were retained for the detailed design. However, the rotor blade 
number and rotor configuration were reevaluated during the detailed design 
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of the rotor using computed boundary-layer characteristics As a result 
of a more detailed evaluation of alternatives, a tandem blade row con- 
figuration with 12 blades per row was selected based on a further im- 
provement In rotor efficiency. Table II presents a comparison of the 
previously tested design and the selected redesign in terms of various 
efficiency expressions and power output. Table III is a reproduction of 
the computer output of the program used to predict the design-point char- 
acteristic of the selected design; the output is considered to be self- 
explanatory. The vector diagrams for the row inlet or exit design sta- 
tions are presented in Figures 7 and 8. 

The 11 initially identified independent variables were reduced 
to 9 by the introduction of two ratios, and were 

treated as fixed values. Subsequently, the ratio was varied for 

the near-final design standard. Each of the variables is briefly dis- 
cussed below using the design of Reference i as the reference turbine in 
the comparisons. 

1. Nozzle Row Inlet Diameter, Compared with the existing 

turbine, this geometric variable was increased by 15 per 
cent. While the selected prediction method indicated still 
further increases in efficiency with diameter, the rate of 
increase was relatively small at the selected value. The 
increasing frictional loss was considered sufficiently large 
that the choice of an optimum solidity for the nozzle row 
would become more critical. Hence, the nozzle inlet diame- 


ter was not increased to the value indicated by the stage 
optimization which implicitly assumes that an optimum so- 
lidity could be selected. 

2. Rotor Inlet Diameter, The diameter was increased by 

6 per cent. The aerodynamic optimum, that is maximum over- 
all tota 1 -to-stat i c efficiency, was predicted to occur at the 
6 per cent greater value °f^g- The diameter increase im- 
plies the use of backslope blading; the actual value of 
rotor blade inlet angle was to be obtained during the de- 
tailed design phase. The stress level increases associated 
with the change were acceptable. 
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3. Nozzle Row Inlet Axial Width, ( 5 ^ . This width was increased 
by 45 per cent. The increase leads to a greater increase 

in velocity across the row. Although the prediction proce- 
dures indicated still further small increases with increase 
of width, the final value was selected after a review of 
possible row inlet blockage effects which could be expected 
to be more significant as the inlet flow angle increased 
still further beyond 55 degrees, 

4. Nozzle Row Inlet and Rotor Row Inlet Axial Widths, a nd 

. These widths, which were assumed to be equal in the 
preliminary design analysis were increased by 4 per cent. 

The new values were predicted to be optimum with increases 
in nozzle row losses offsetting decreases or rotor losses 
with any further increase in the flow path width. 

5 . Rotor Row Exit Outer Diameter, 33-^. This value Is unchanged 

5t 

from the reference turbine, since it corresponds to a near 
optimum with rotor and diffuser losses changing at ap- 
proximately equal but opposite rates in the vicinity of the 
unchanged value of"^^. 

6. Rotor Row Exit Inner D i ameter , This value was also 

unchanged. When investigated keeping the rotor exit annulus 
area constant by also varying the outer diameter, the origi- 
nal value was concluded to be near optimum. 

7. Number of Nozzle Row Blades, jZ . The selected number of 
nozzle blades was 15 with this number predicted to corres- 
pond to an optimum solidity, (In the final evaluation of 
the efficiency of the redesign, one of the alternative pre- 
diction procedures indicated that this number of blades may 
correspond to a solidity in excess of an optimum.) 

8. Number of Rotor Blades, . The initial recommendation 

of the NREC investigators was for 10 full blades and 10 
splitters. However, the NASA personnel's evaluation of the 
selection resulted in a change to a tandem row configuration 
which was also approved by NREC on the basis that the turbine 
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9. 


was to be rig tested and a tandem row configuration might 
result in advances in turbine technology. While the pre- 
diction procedures had indicated a small reduction in rotor 
losses for a tandem row configuration with an appropriate 
number of blades per row, the lack of experimental data to 
confirm the choice of such a configuration was regarded as 
a relatively high risk approach; performance data for a tandem 
row rotor, however, were agreed to be of considerable inter- 


est. 

Rotor Axial Length, The rotor axial length was reduced 

by 12 per cent. While not shown to be particularly signifi- 
cant in terms of rotor performance in the preliminary analy- 
sis, any further reduction of length could have resulted 
?n increased blade lean in the vicinity of the maximum cen- 
trifugal root stress. 

Vaneless Space Ratio, This value was unchanged at 

1.03* While a reduction was predicted to be marginally 
beneficial to stage performance as a result of vaneless 
space losses, the ratio was maintained at the value selected 
for the reference turbine to avoid possible effects asso- 
ciated with nozzle and rotor interaction with any closer 


spacing of the nozzle row and the rotor. 
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AERODYNAMIC DESIGN 

Rotor Design 

Objectives and Approach 

The principal objective of the aerodynamic design is to define 
the meridional flow path and blading geometries of a tandem row rotor.. 
Based on an evaluation of the results of flow analyses, the geometric 
specification of the rotor must, of course, be consistent with the over- 
all design intent of designing a high performance stage. 

The preliminary design analysis had established the principal 
dimensions of the rotor together with the inlet and exit flow conditions* 
The detailed design phase consisted of an iterative procedure making use 
of the computerized design methods outlined in Appendix I. Although the 
preliminary design analysis had led to the selection of 12 blades in each 
row, during the detailed design it was established that either more blades 
or a greater axial length would be required* Of the two alternatives, the 
former was selected and the number of blades increased to 16 in each row. 
The basic reason for the change was that the loading over the rearward 
portion of the first row had to be reduced below that originally targetted 
in order to avoid flow separation from the suction surface of the first 
blade row; any substantial separation might be expected to offset any 
other advantage associated with a tandem row design. The loss in load- 
ing associated with satisfying the blade end condition of equal velocities 
on both suction and pressure surfaces for the first row while maintaining 
an acceptable diffusion rate on the suction surface had been underesti- 
mated in the preliminary design analysis comparison of tandem row and 
full blade configurations. 

For the revised number of blades, the design objective in terms 
of the surface velocity distribution is shown in Figure 9 ; the velocities 
shown are for an outer casing streamtube. The finally selected meridional 
flow path through the rotor was based initially on an analysis in which 12 
blades were assumed. The objective loading for 12 blades is not illus- 
trated herein. However, to a good approximation it was that which would 
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be obtained by loading the first blade row so that the composite loading 
would approach that of a single blade. The rapid unloading near the 
trailing edge of the first row was subsequently determined to be unachiev- 
able without a massive separation of the flow. Since the mean relative 
and absolute velocities for the design would not be qreatly influenced by 
the change in the number of blades, the meridional flow path geometry was 
not redefined following the increase from 12 to 16 blades per row. 

Rotor Meridional Flow Path 

The selected rotor flow path is shown in Figure 10. The rotor 
tip diameter, the Inlet axial width, the annulus diameters at rotor exit, 
and the axial length of the rotor are as determined by the preliminary 
analysis. The figure shows the axial location of the split between first 
and second rows. The hub line of the rotor is a spline curve which matches 
the slope of stator wall at rotor inlet. Figure 11 shows a comparison be- 
tween the redesigned and the original stages. It will be seen that only 
the rotor exit annulus dimensions are unchanged and that the major rotor 
flow path changes are in the diameter of the rotor, axial length, and in 
the details of the inner and outer contours. It can be seen that the 
redesigned flow path has reduced curvatures along the inner and outer con- 
tours. 

The specifications for the rotor flow path are given in Table 
IV. This table contains the r, z coordinates of the hub and the outer con- 
tour of the rotor; actual stationary casing dimensions are obtained by the 
addition of clearance to the rotor. The pairs of hub and outer contour 
points given in the table correspond to the ends of straight-line elements 
of the blading. 


Fi rst Blade Row 

While the blade camber line definition for a streamline near 
the outer casing could be computed from the meridional velocity data of 
preceding flow analyses and the targetted loading distribution, the load- 
ing along the hub streamline was constrained by mechanical design consid- 
erations. The orientations of the stra i ght-1 i ne bl ade elements defining the 
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blade were constrained to produce a mechanically acceptable blading stan- 
dard; the first row trailing edge line is approximately radial. As a re- 
sult of the constraint, the split of loading between the rows along the 
hub line was initially a variable in the design iteration. The major 
portion of the flow path was derived using the ANALS computer program as 
the design tool. Blade inlet and exit regions, however, were finalized 
using Program BL2BL as the flow analysis tool. 

The finally selected blade geometry is defined in Table V. This 
table is reproduced from computer output; points of corresponding number 
define straight-line elements through the mean camber line of inner and 
outer streamlines with the thickness distribution defined by cones having 
the specified thickness at the inner and outer contour end points. Blade 
inlet and exit angles are implicit in the specification of the blade camber. 
The calculated blade inlet angle varies from approximately -10 degrees at 
the hub to -13 degrees at the outer contour; thus, the amount of backs lope 
in the rotor inlet is relatively small with the rotor operating in the range 
of -20 to -17 degrees of incidence at the design point. The calculated 
deviations at first row exit are relatively small as a result of the specifi- 
cation of relatively lightly loaded trailing edge regions. 

The calculated velocity distributions for two streamlines near 
the outer contour and near the hub are shown in Figures 12 and 13. These 
plots of blade surface velocities as a function of meridional distance 
combine data from Programs BL2BL and ANALS with the former used for the 
inlet and exit regions and the latter for the major portion of the flow 
field where the assumptions related to a channeled flow are considered 
valid. The figures demonstrate that diffusion of suction surface ve- 
locities has been controlled and that pressure surface diffusions have 
been restricted to the immediate vicinity of the leading edge. 

Second (or Tandem) Blade Row 

From Figure 10 it can be seen that the second row can be con- 
sidered as axial blading to a good approximation. Although row inlet and 
exit annulus dimensions correspond to their being part of a radial inflow 
turbine, the radial components of velocity at both stations are relatively 
minor compared to the axial and tangential components of velocity. Hence, 
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if a satisfactory distribution of loading and surface velocity distribu- 
tion could be obtained ignoring the radial components, it is reasonable 
to expect an acceptable distribution when the resultant geometry was sub- 
jected to standard three-dimensional flow analysis. 

The program used to obtain the geometric specifications of two 
constant radius sections (one near the casing and one near the hub) is an 
NREC computer program based on the Stanitz prescribed velocity distribu* 
tion method of Reference 10. To select the prescribed distributions of 
the datum section, the following guide lines were established: 

1. The inlet flow angles were those corresponding to the axial 
and tangential components of velocity at the first row 
exit at the appropriate radii. 

2. The exit angles and Mach numbers were as previously estab- 
lished from the preliminary design results. 

3. The velocity distributions were basically those previously 
established as objectives. 

4. The leading edge loading would be such that there would be 
a negligible effect on the upstream flow angles which would 
be established by the first row. 

5. The loading distribution over the rearward portion of the 
aerofoils would be such as to avoid any separation. Hence, 
the deviation of the flow would be predictable by potential 
flow methods. 

6. The stack of the sections together with the area taper ratio 
should be consistent with the over-all mechanical require- 
ments. 

Profiles meeting the above requirements were obtained; the profiles to- 
gether with the prescribed Mach number distributions used to generate 
them are shown in Figures 14 and 15- Using these basic profiles and 
assuming straight-line generation of the complete second blade row, the 
flow field analysis was completed using Programs ANALS and BL2BL. The 
calculated blade surface velocity distributions for two streamtubes near 
the casing and near the hub are shown in Figures 16 and 17; these velocity 
distributions indicate that the design intent was achieved. 
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The geometric specification of the second blade row is presented 
in Table VI. The geometric data are given for two streamlines with cor- 
responding camber line points defining straight-line elements and normal 
thickness defining conical tapers. 

In the table of angular coordinates, the reference angle for the 
second row is arbitrary; for the final design specification contained in 
the rotor drawing the leading edge of the second row is located approxi- 
mately midway between the trailing edges of the first row. 

Although the two blade rows were designed independently in terms 
of the blade-to-blade flow field calculations, the design approach to the 
rearward portion of the first row and the leading edge region of the second 
row was such that there is reason to believe the calculated surface ve- 


locity distributions are acceptably accurate. The blading, when viewed 

in terms of the streamline blade angles, & = £an^/e4 s shows acceptable 

' avr) 

continuity between the rows; this can be interpreted as having designed 
blading with small flow deviation from the first row and small incidence 
at the second row. Figure 18 presents a blade angle plot for both rows 
for a streamline near the outer casing. 


Nozzle Desi 


The detailed definition of nozzle vanes and meridional flow 
path geometries and the calculation of free stream velocity distributions 
were accomplished using the approach described in Appendix I. An initial 
vane geometry design was generated from a prescribed distribution of load- 
ing. Using this same specification and the meridional flow path shown 
in figure 19 as input to the flow analysis program ANALS, acceptable dis- 
tributions of surface velocities were obtained. 

The surface velocity distributions shown in Figure 20 for the 
nozzle vanes were obtained using the more rigorous method of Program 
BL2BL. Both suction and pressure surfaces show continuous acceleration 
of the flow from inlet to exit. The calculated deviation at the trailing 
edge is 0.4 degrees, and the average flow angle at the leading edge of 
the rotor was calculated to exactly equal the required value of 74.1 
degrees . 
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The selected nozzle geometry is illustrated in Figure 21; the 
profile consists of a 6 per cent thick NACA-63 airfoil superimposed on a 
camber line computed from the preselected vane loading distribution- The 
geometric specification of the vanes is given In Table VII; the vane blad- 
ing is of constant section in the axial direction. 

Exhaust Diffuser Design 

Since the dimensions of the turbine rotor exit annulus are un- 
changed (the exhaust diffuser was constrained to match the testing de- 
livery system) and the total length of the existing diffuser was considered 
adequately near optimum, the design of the exhaust diffuser consisted pri- 
marily of the calculation of flow area as a function of distance from the 
rotor. 

NASA test data of Reference 11 have shown a diffuser designed 
for a linear static pressure rise to improve the performance of the exist- 
ing BRU stage. Hence, the diffuser of the subject turbine has also been 
designed for a linear static pressure rise over the major portion of its 
length. A constant diameter was selected for the diffuser center body 
since any residual swirl will adversely affect the performance of an annu- 
lus diffuser having a reducing value of inside diameter. Hence, the dif- 
fuser aerodynamic design effort consisted of the calculation of the outside 
diameter for assumed values of total and static pressure as a function 
of axial length. Since the available length exceeded that necessary for 
the generally accepted optimum diffusion angle, an initial section of 1.25 
in (3*175 cm) was designed for a constant static pressure, and the final 
exit section of 1.0 in (2.5^ cm) was selected to provide a smooth transi- 
tion to the existing delivery ducting. Over the linear static pressure 
portion of the diffuser the static pressure was assumed to vary between 
the diffuser inlet and exit static pressure values determined during the 
preliminary design analysis (see Table III). Similarly, the linear drop 
of total pressure through the diffuser is between the total pressure val- 
ues established during the preliminary design. 

The geometric specification of the diffuser is given in Table 

VIII. 
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Inlet Scroll Des I gn 

The scroll was designed to match an existing 5 in diameter in- 
let duct. The design procedure used is that detailed in Appendix i with 
the type of scroll geometry illustrated in Figure 22. Since the nozzle 
row inlet dimensions have been changed considerably in the redesign tur- 
bine compared with those of the reference turbine, a major change in the 
scroll detail was required. It was decided to design for a more compact 
scroll. As a result of the more compact design the tangential distance 
between the center line of the inlet ducting and the axis of the turbine 
was reduced. This change can be readily accommodated on the turbine test 
stand. 

Even though the ssure loss in the inlet scroll is predicted 
to have little effect on tnc over-all performance of the stage, the pre- 
viously computed loss was included in the calculation of scroll areas. 

The area and mean radius schedules for the scroll were calculated to be 
consistent with the design intent of a circumferentially uniform distri- 
bution of static pressure and radial component of velocity at nozzle row 
inlet. 

The geometric specification of the inlet scroll is given in 

Table IX. 

Boundary-Layer Analysis 

Two-dimensional boundary-layer calculations were undertaken for 
representative surfaces of the rotor and the nozzle row. The procedure 
used was that outlined in Appendix I; the computer program used was Pro- 
gram BIANCA of Reference 12. For the rotor, eight surface elements were 
analyzed corresponding to pressure and suction surfaces of hub and casing 
streamtubes for both blade rows. For the nozzle row, two stream surfaces 
correspondi ng to the pressure and suction surface were analyzed. 

These bounda ry- 1 ayer calculations ignore the effect of cross 
flows; that is, the migration of low momentum boundary layers in directions 
other than the free stream direction of the inviscld flow field. However, 
two secondary flow parameters, Cp and 2f , were also calculated to 
provide a qualitative indication of the magnitude of secondary flow 
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effects. The first of these parameters serves an an index for end wall flow 
and, in the case of the rotor, tip leakage. The second serves as an Index 
for boundary-layer flow migration on blade surfaces. 

R otor Data 

Figures 23, 24, 25, and 26 present the results of the boundary- 
layer analysis for the first row of rotor blading. Each figure contains 
two sets of normalized momentum thickness and shape factor together with 
the relevant secondary flow parameter for the flow region considered. 

Figures 27, 28, 29, and 30 present corresponding data for the second 
row. In only one case was separation predicted, this being for the suc- 
tion surface of the first blade row at the hub streamtube. However, since 
the secondary flow will be away from this region of the blade, separation 
may not in fact occur. Conversely, the flow of low momentum fluid to ad- 
jacent streamtubes could cause the separation region to extend. In either 
event flow separation is not likely to result in any large efficiency 
penalty since the rearward portion of the first row was intentionally 
1 i ghtly loaded. 

While the data presented in Figures 23 through 30 serve prin- 
cipally to illustrate the behavior of the boundary layers In the rotor, 
the data have been used for the design-point performance prediction dis- 
cussed later. One tentative conclusion which can be drawn from the 
boundary- layer analysis is that the rotor design still does not consti- 
tute an optimum configuration even though it is expected to outperform 
the reference rotor. Based on a comparison of data from the first and 
second rows, it would appear probable that a configuration using more 
blades in the first row and fewer blades in the second row would result 
in lower rotor losses. Thus, a tandem row design in which the number of 
blades in the first row is twice that of the second row is recommended 
for any further investigation of tandem row radial inflow turbine rotors. 

Nozzle Row Data 

Figures 31, 32, and 33 present results of boundary-layer analy- 
ses for the nozzle row. In the case of the nozzle row the data are expected 
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to be approximately equal for the forward and rearward end walls. As was 
to be expected from the previously discussed surface velocity distribu- 
tions, relatively low values of normalized momentum thickness and shape 
factor are predicted for the surfaces of the nozzle row. However, the 
data have been used to calculate a row total pressure loss from which it 
is inferred that the selected solidity of the nozzle row is in excess of 
an optimum value. 


Design-Point Efficiency Estimate 

Following the completion of the aerodynamic design of the noz- 
zle and rotor rows and the boundary- 1 ayer analyses discussed above, a 
reevaluation of design-point performance is possible. Total pressure 
loss coefficients for the nozzle row and the rotor were obtained using 
the procedure presented in Reference 13* in turn, these losses can be 
expressed as efficiency decrements for comparison with the values obtained 
at the completion of the preliminary design phase. The following tabula- 
tion presents the comparison and also includes the diffuser decrements 
which were not reevaluated. 


Efficiency Decrement 


Component 


Prel i mi nary Des i gn 


Final Design 


Stator 
Rotor 
Diffuser 
Tota 1 


0.0161 

0.0402 

0.0137 

0.0700 


0.0341 

0.0440 

0.0137 

0.0918 


The stator comprises the inlet scroll, the nozzle row, and the vaneless 
space between the nozzle row exit and the rotor inlet. The above total 
decrements yield efficiencies based on over-all scroll inlet total pres- 
sure to diffuser exit static pressure of 93 end 90.82 per cent, respec- 
tively. If it is assumed that the more rigorous loss analysis undertaken 
following the completion of the detailed aerodynamic design is a more re- 
liable estimate of loss than that undertaken during the final design, the 
performance improvement compared with the existing design will be 2 rather 
than 4 percentage points. While the higher loss in the rotor can be traced 
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to the increase in the number of rotor blades, the more significant change 
is that in the prediction of the stator system (and in particular the noz- 
zle row) loss and Is less readily explained. Which of the two alternative 
prediction procedures will predict most accurately the total pressure loss 
in the nozzle row will have to be determined from turbine rig tests. It 
should be noted that both prediction methods are the result of adapting 
existing methods for other turbomachinery components to the nozzle row 
of a radial inflow turbine and the particular case of a row across which 
there Is a significant reduction in blade height. While the prediction 
procedure used in the preliminary design indicated that 15 blades was an 
optimum for the nozzle row, the loss derived from boundary- 1 ayer analysis 
can be interpreted as an indication that this number exceeds the optimum. 
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MECHANICAL DESIGN 


I ntroduct I on 

The detailed mechanical design was carried out such that the 
research parts are consistent with the existing bearing housing and stand 
and will meet the package and interface requirements of the BRU, as de- 
fined in Reference 2. The various procedures adopted in the mechanical 
design analyses are documented in Appendix II. 

The turbine research parts were designed for ease of assembly 
so that alteration and/or replacement of the rotor, stator blade row, and 
exit diffuser are facilitated. An outline drawing, showing interface lo- 
cations and useful over-all dimensions, is provided in reduced size form 
In Appendix III, which also contains reproductions of all drawings referred 
to below. The stress analysis of the rotor was performed to establish an 
adequate overspeed safety margin under hot design conditions. 

A critical speed and bearing load analysis of the rotor system 
was performed to insure that critical speeds are out of the 10,000 to 
35,000 rpm range for performance testing and that the bend shaft criti- 
cal occurs above 120 per cent of the maximum operating speed. The burst 
speed of the impeller disk at design temperature was verified to be above 
140 per cent of design speed. Peak steady-state dynamic rotor stresses 
at 120 per cent of design speed were calculated to be below the 0.2 per 
cent yield strength for the rotor disk material, including the equiva- 
lent stresses in the blade attachment areas. Extreme operating condi- 
tions were evaluated to assure that blade tip clearance, stress, and 
balance conditions are acceptable. 

Blade stresses due to rotation, aerodynamic loading, and tem- 
perature were evaluated and a Goodman diagram plotted to determine the 
margin of safety for alternating stress. A Campbell diagram was plotted 
to confirm nonresonant operation from 40 to 120 per cent of design speed. 

The design was provided with all instrumentation required with 
locations consistent with the specified requirements. 
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Turbine Scroll and Nozzle Design 

The turbine scroll and nozzle ring were designed to meet the 
test operating conditions specified for component testing which included 
turbine inlet temperature up to 150 deg F (338 deg K) and an inlet pres- 
sure range from 10 microns of mercury to 60 psia (414 kN/m ). The inlet 
flange was designed to interface with the stub-end flange as specified in 
AiResearch Drawing No. 699770. Calculations of both maximum stress and 
axial deflection were performed, and results indicated values to be well 
within acceptable limits for the materials selected, 

Hastelloy X was selected as the material and a welded sheet 
metal design was established in order to minimize tooling changes and 
maintain compatibility with requirements for the BRU as defined in Ref- 
erence 2, Minor changes in the inlet line location were coordinated 
with the NASA Project Manager in order to incorporate a scroll design 
which minimizes the vaneless annular space upstream of the nozzle vane 
inlet and hence ha^ a smaller maximum diameter than that shown in Refer- 
ence 2, The revised location of the scroll inlet flange is defined in 
the research turbine assembly outline (Drawing Nc, 1 1 59-R-003) . 

The turbine scroll inlet section between the inlet flange and 
the scroll entrance station was carefully designed to provide a smooth 
area transition from the 5 in (12,7 cm) diameter circular inlet area to 
the initial typical cross-sectional area at the 0 degree scroll entrance 
station as shown in Figure 34. Details of the transition can be seen in 
the scroll fabrication assembly drawing (Drawing No. 1 1 59-R-006) . 

Primary emphasis was placed on ease of assembly and alteration 
or replacement of component parts as can be seen in the cross-section 
provided in the research turbine package assembly drawing (Drawing 
No. 1 1 59-R-0Q2) . The scroll is bolted directly to the bearing housing 
with a shim ground to size and placed under the scroll mounting flange 
to establish the rotor operating clearance of 0.009 to 0.011 in (0.0229 
to 0.0279 cm). Provisions were made such that the turbine rotor, turbine 
nozzle ring, and turbine exhaust stationary shroud, including the diffuser 
section can all be installed or removed without disturbing the scroll as- 
sembly. 
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Rotor Mechanical Design 


Rotor Steady-State Heat Transfer Coefficients 


Rotor steady-state heat transfer coefficients, Biot number, and 
adiabatic wall temperatures are required to calculate rotor metal surface 
temperatures which are used, In turn, to predict the temperature distribu- 
ting within the rotor. For analysis purposes, the rotor is subdivided 
into discrete surface areas, as indicated in Figure 35, and heat transfer 
coefficients are calculated based upon known and/or assumed boundary con- 
ditions. For the disk surfaces, such as BC, CD, and EF in Figure 35, the 
heat transfer coefficient depends upon the Reynolds and Prandtl number 
and is calculated by means of equations derived empirically for high ve- 
locity flow, the general expression for which is as follows: 
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For the areas where blades exist, such as surface DE in Figure 
35) the surface is treated as a finned surface and the additional heat 
transfer is accounted for by a correction to the standard heat transfer 
coefficient based upon the appropriate area ratio as indicated in the 
above equations. 

In the rotor attachment area, AB in Figure 35* the heat trans- 
fer coefficient is assumed to be simply the ratio of metal conductivity, 
K Me L> divided by the length of the attachment surface or heat transfer 
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path, L » and the heat flux is calculated in this area from the boundary 
temperature conditions, one of which is assumed to be the temperature 
of the shaft at the bearing location. 

After the heat transfer coefficient has been computed, the heat 
transfer between the fluid and the solid body is further described in the 
appropriate areas of Figure 35 by the standard dimensionless grouping of 
variables further identified as the Biot number and defined as 



( 6 ) 


where ^g^is effective heat transfer coefficient 
Y^p is outside radius of disk 
and K^ a ^is metal conductivity. 

Finally, adiabatic wall temperature is reevaluated from the 
genera] equations 


■aw oo ' / T 

■fr = (P ,-/ 2 


(7) 


The adiabatic wall temperature , , is based upon the temperature as- 
sumed by an insulated plate when exposed to a high velocity hot gas stream. 
The recovery factor, » assumes a turbulent boundary layer associated 

with high velocity flow and is based upon Prandtl numbers in the range 
between 0.5 and 2.0. The variables in the above equation are local static 
temperature , , local relative velocity, W , fluid speci f i c heat at 

constant pressure, Qp , and Prandtl number, Py" . 

For the blades portion of the disk, DE, relative velocity and 
static temperatures are evaluated at points along the mid-streamline of 
the blade passage. Relative velocity on the back, or unbladed portion of 
the disk, BC, is assumed to be one half of the local disk tangential ve- 
locity. 
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A summary of the results of these calculations is provided in 
Table X, which lists the Biot number and adiabatic wall temperature as a 
function of radius for the various surface areas of the turbine rotor. 

Calculation of Temperature 
Distribution in Rotor 

The rotor adiabatic wall temperature and Biot number data were 
used to determine the temperature distribution both on the surface and 
internal temperatures of the rotor. The results of this analysis are 
plotted in Figure 36, which shows constant temperature lines plotted on 
a cross-section of the rotor disk. The analysis was performed using 
Program THERMAL STRESS of Reference 4. 

Thermal and Mechanical Disk Stress 

An analysis of disk stresses and deflections was also performed 
with the use of Program THERMAL STRESS, which produces results based upon 
the geometry, operating speed, and temperature conditions provided. Both 
radial and tangential stresses were calculated and these are shown as a 
function of radius in Figures 37> 38, and 39 for 100 per cent, 120 per 
cent, and 140 per cent speed, respectively. Radial and tangential stress 
distributions are also shown in Figures 40 and 4l , respect i vel y ,. for the 
100 per cent speed condition. All of the above information is provided at 
design operating temperatures of the BRU rather than the turbine stage re- 
search package. 

The selected material for the rotor is Inconel Alloy 713LC; 
this was the material selection for the original turbine also. Before 
making the decision to retain the material specification of Reference 2, 
other materials were also evaluated. Table XI lists the properties of 
the materials considered. Reference to the properties of 713 LC shown in 
Table XI and the stress data in Figures 37 through 41 will show that 
there are no advantages associated with a change in the rotor material 
nor is a change necessary to meet the turbine life objective. It is evi- 
dent that the rotor design easily satisfies the requirements that the 
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maximum steady-state dynamic design stress at 120 per cent be less than 
the 0.2 per cent yield strength and that the burst speed of the rotor 
testing, based upon material ultimate strength, shall be above 140 per 
cent of the design speed. Equivalent stress, as usually defined, is 
shown in Figure 42 as a function of radius for both the front and back 
faces of the disk at 120 per cent speed and design operating temperature. 

Calculated values of average tangential stress are provided in 
Table XII for 100, 120, and 140 per cent operating speeds for design op- 
erating temperature conditions. Also shown in the table are calculated 
values for average tangential stress at room temperature for a disk of 
alloy 7 1 3 LC material at 120 per cent speed and aluminum alloy 7079-T6 at 
100 per cent speed. Axial and radial disk tip deflections calculated 
for the five cases abo'" are also included in Table XII. 

Axial and radial deflection, as a function of radius, are pro- 
vided in Figures 43 and 44, respectively, for the 120 per cent speed con- 
dition at operating temperature. These figures also include data for 
equivalent deflections at the lower level to be used in turbine rig tests. 

Based on cost and delivery estimates assuming three or less 
rotors would be procured, a machined rotor is recommended. 

Blade Root Stresses 

Blade root stresses due to centrifugal loading were calculated 
using Program BLDR0T of Reference 14. The blade geometry is considered 
as a series of beams of rectangular cross section, pinned at the shroud 
line to adjacent beams but otherwise unrestrained except at hub; rela- 
tively conservative values of computed stress are expected from this ap- 
proach. Blade thickness distributions used at the hub and shroud of the 
first and second blade rows, respectively, are shown in Figures 45 and 46. 
Calculated blade root bending and tensile stresses at 120 per cent speed 
are presented in Figures 47 and 48 for the first and second blade rows, 
respect i ve ly . 
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Blade Natural Frequency 

Blade leading and trailing edge frequencies were calculated 
based upon plate theory. The plate model was assumed clamped on two sides, 
representing its attachment to the rotor hub and to the adjacent blade sec- 
tion or panel. Results are summarized In Table XIII; the lowest frequency 
was, as expected, calculated to be that for the second blade row trailing 
edge, A Campbell or interference diagram is plotted in Figure 49; the blade 
natural frequencies are shown relative to multiples of the rotor “once-per- 
revolution" frequency, in cpSn 

Blade Vibratory Stress Levels 

For the purpose of calculating blade vibratory stress, blade 
aerodynamic loading was calculated from the blade static pressures at hub 
and shroud and assuming a linear pressure distribution between these two 
points. An amplification factor based on the method of Trumpler and Owens 
(Ref 15) was used, with an appropriate stress concentration factor and the 
blade aerodynamic bending stress, to obtain a design value for the blade 
vibratory stress. The calculated results are shown in Table XI II and on 
the Goodman diagram provided as Figure 50. While predicted vibratory 
stress levels in the second row are not insignificant, they are never- 
theless acceptably low. 

Flexural Vibrations 

A shaft flexural vibration analysis was performed on the rotor 
system illustrated in Figure 5 1 -- The system consists of the turbine, 
the rotor shaft mounted in two bearings, and half the aluminum disk coup- 
ling mounted on the output end of the shaft. The critical speed analysis 
was performed using an NREC computer program, Program HOLZER, which is 
based on standard procedures presented in Reference 16 and outlined in Ap- 
pendix M. The analysis was completed for a range of turbine end bearing 

c 

stiffness while maintaining a rear bearing stiffness of 2.5 x 10 1 bf per 

in (4,38 x 10' N/cm) (as assumed in Ref 2). The results of the calcula- 
tions are shown in Figure 52 which illustrates that there are no critical 
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speeds in the operating range from 10,000 to 43,200 rpm for an assumed tur- 

4 4 

bine end bearing stiffness of 3 x 10 Ibf per in (5»25 x 10 N/cm) (the 
stiffness assumed in Ref 2). 

An analysis was also performed of a similar rotor system but as- 
suming an aluminum turbine disk. As can be seen in Figure 52, the first 
rotor critical speed now occurs at 14,000 rpm, which is within the speci- 
fied operating range and hence unsatisfactory. An aluminum alloy rotor 
was considered for the research package because stress and temperature 
levels would permit its use, reduce manufacturing costs, and facilitate 
manufacture. 

Shaft Attachment Methods 

The rotor will be attached to the shaft by means of a heavy 
shrink fit, as indicated in turbine wheel assembly Drawing No. 11 59-E-005. 
Based on considerable experimental evidence, the heavy interference fit 
has more than adequate torque capacity and axial holding stability for 
the research application planned and will be used without radial pins of 
any kind to avoid the installation and operating problems which may ac- 
company such devices. The stress concentrations resulting from radial 
holes in the highly stressed attachment area are also avoided. 

Beari ng Analysi s 

The adequacy of the existing bearing of the turbine research 
package was investigated. It was concluded that the bearing will meet the 
design objective of a time between overhaul (TB0) of 300 hours minimum. 

This conclusion is based on the calculated value of the expected minimum 
B ]0 life of the bearing. 

The bearing used in the turbine research package is described ; 
in Reference 2. The description consists of the major characteristics of i 
the bearing and a source control drawing. Simply stated, it is a Type 
7104 angular contact ball bearing. As such, its B^q life should be ob- 
tained using the method established by the Anti -Friction Bearing Manu- 
facturers Association (AFBMA) . This method requires the calculation of 
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the thrust load and the radial load acting on the bearing. The details 
of the calculations are given below. 


( 8 ) 


where 


Pf 

Pb 

Vrv 


and 


Thrust Calculation 

The calculation of the net thrust on the wheel is based on con- 
siderations of the wheel as an open system. The axial forces acting on 
the wheel arise from the frontal and back pressures of the fluid and the 
change in axial momentum of the fjuid. Hence, 

Thrust = f \pj_ 2Ti'ro!r'+ rnc^ 

° J % 

is the frontal pressure 

is the back pressure of the fluid 
is the tip radius of the wheel 
i s the mass f low 

Cy^^is the axial velocity at the turbine exit, 

A net thrust in the direction of the exit flow of 74»9 1 bf 
(333.2 N) was obtained for the fluid properties and the design flow con- 
ditions. More specifically, a mass flow of 0.748 1 bm per sec (0.3395 kg/ 
sec), an exit axial velocity of 24 1 fps (73°5 m/sec), and the pressures 
shown in Figure 53 were used to evaluate the thrust. The frontal pressure 
in Figure 53 corresponds to static pressure at the shroud for 1.730 in iSVig 
2.634 in (4.394 6.690 cm), a linear variation between the hub 

and shroud values of the exit states pressure for 0=911 in 1.730 in 

(2.3i4^y^ 4.394 cm), and the exit static pressure at the hub for 0 
0.911 in (2,314 cm). The back pressure, on the other hand, corresponds 
to a variation from the rotor inlet static pressure defined by 

k . Ekh 

w- %r 

where^^ is the fluid density and W is the relative velocity. Assuming 
that the fluid in back of the wheel rotates at one half of the wheel speed, 

CO , the pressure variation may be rewritten as 

dp I Po»V 

"Sr 4 So 10 

The density is evaluated at the turbine inlet conditions. 
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Radial Load Calculation 

The radial load exerted on the bearing is obtained using the 
n 2g" approximation. That is, the reaction of the bearing to the dynamic 
radial load is assumed to be less than or equal to the reaction to twice 
the static load. The reaction of the thrust bearing, can be ob- 

tained from a summation of the forces and moments. 

A reaction of 11.1 lbf (49.4 N) was obtained based on a weight 
of the turbine wheel, w r , of 3-868 lb (1.75 kg), a distance between the 
centers of gravity of the wheel, and the thrust bearing of 2,371 in (6.022 
cm), and a distance between the centers of gravity of the bearings of 
4.006 in (10.175 cm). Distributing twice the weight of the shaft evenly 
between the bearings results in additional reaction of 2.1 lbf (9*34 N) 
or a total radial load of 13.2 lbf (58.72 N) for the thrust bearing. 

Expected Mi nimim q Life 

The AFBMA method of calculating B|q life states that 

3 

minimum life in hours = 500 x (service factor) 


where 


service factor 


rated capac i ty 
equivalent load 


The rated capacity of a Type 7104 bearing at the required speed of 36,000 
rpm is obtained according to AFMBA procedures from 


rating at required speed = rating at 1000 rpm x 3 


1000 rpm 
requ i red speed 


Hence, a rated load at 38,000 rpm of 2 1 4 lb (952 N) is obtained based on 
a rating at 1000 rpm of 710 lb (3158 N) . The equivalent load of any 
7100 series bearing is obtained by selecting the maximum value from 

£=• 0 . 51 VR + T or E = 1.31 VR 

where EL is the equivalent load 

R is the radial load 

T" is the thrust load 

and V '• 1.0 for an inner rotating race. 
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Using the results previously obtained and the preferred bear- 
ing preload of 60 lbf (267 N) , an equivalent load of 142 lbf (632 N) is 
obtained. Hence, from calculated rated capacity and equivalent load, the 
minimum life is I 67 O hrs. 

Diffuser Design 

The diffuser assembly is shown in diffuser machining Drawing No, 
1159-D-010. The Installation of the diffuser is shown in both the tur- 
bine research package assembly Drawing No., 1 1 59-R -002 and the turbine re- 
search package outline Drawing No. 1I59“R"003. Instrumentation pads 
for the total pressure traverse at inlet and exit are provided in addi- 
tion to static pressure taps on the outer wall and on the center body at 
both inlet and exit. The center body is supported axially and radially 
and maintained cov-c.in.ric with the outer contour by two sets of three 
aerodyne,, ii ca 1 1 y shaped struts installed near the inlet and exit of the 
<ii rfuser, as shown on the referenced drawings. 

It is assumed that an external support exists and will be pro- 
vided and used in the final assembly to prevent distortion of the scroll 
and nozzle assembly, due to the large overhung mass of the diffuser. The 
exhaust adapter Drawing No. 1159**D-011 was designed to interface with the 
existing laboratory exhaust system in both length and diameter, as speci- 
fied. 

Final Assembly and Installation 

Copies of all detail drawings produced, plus the assembly, out- 
line, and layout drawings are provided in Appendix III. The assembly draw- 
i ng includes a complete parts list, also shown separately, with miscellane- 
ous tools required In Table XIV. Assembly procedures are covered in the 
notes on appropriate drawings and all instrumentation locations are estab- 
lished in either the detail drawings, the assembly, or the outline drawing. 
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CONCLUSIONS AND RECOMMFNDATI ONS 


Cone 1 us i ons 

1. A complete radial inflow turbine stage comprising an inlet 
scroll, a vaned nozzle row, a rotor, and an exhaust diffuser 
has been designed and drawings for all new parts prepared. 

The stage is a redesign for a BRU turbine previously tested 
in a turbine research package. The new design is the result 
of an aerodynamic optimization; although the rotor diameter 
has been increased and nonradial blading introduced, the 
increase In stress levels is mechanically acceptable because 
of the relatively low levels of rotor tip speed and the op- 
erating temperature. The rotor mechanical design analysis 
has been completed for the actual BRU operating conditions. 
However, the stationary parts consisting of the scroll, the 
nozzle row, and the diffuser have been designed for the re- 
search package and its lower operating temperatures. 

2. The design-point efficiency of the redesigned stage, based 
on over-all tota 1 -to-sta t i c pressure ratio from scroll in- 
let to diffuser discharge, Is predicted to be between 91 and 
93 per cent representing a 2 to 4 percentage point improve- 
ment compared with the original design. 

3. The uncertainity associated with the predicted stage design- 
point efficiency largely reflects the uncertainty in the 
prediction of the nozzle row loss coefficient. The low loss 
coefficient and high efficiency assessment is based on a 
loss correlation based on axial turbine stator rows which was 
developed for the preliminary design analysis. In this 
method the number of blades is optimized based on a trade- 
off between exit passage aspect ratio and trailing edge block 
age but assumes the row will be designed with near-optimum 
solidity. The efficiency assessment applied to the final de- 
sign standard is based on the calculation of boundary layers 
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with a correction based on wa 1 1-to-blade surface areas. The 
latter efficiency assessment indicates that the solidity of 
the nozzle row is considerably in excess of an optimum. How- 
ever, at this time neither method of predicting nozzle row 
loss coefficients can be regarded as superior to the other, 
since there are no experimental data for nozzle rows of the 
type selected, tnat is, one in which there is a substantial 
contraction of axial width across the nozzle row 

4. The rotor has been designed using a tandem row configuration 
with 16 blades in each of the two rows. The predicted 
efficiency of the rotor is responsible for an approximate 
1.5 per cent increase in stage efficiency compared with the 
original design. One disadvantage associated with the de- 
sign of a two-row configuration highlighted during the final 
design was that a greater number of blades per row would be 
required in a tandem configuration compared with one with 
full rotor blading from inlet to exit. The selection of 
surface velocity distributions to avoid flow separation 
from the suction surface of the first row leads to a limita- 
tion on the blade loading for the first row. 

5- The exhaust diffuser has been designed for a linear rate of 
static pressure rise following an initial section of constant 
static pressure immediately downstream of the rotor. The 
predicted improvement in over-all efficiency of the stage in- 
cludes a 0.9 percentage point improvement for the redesigned 
di ffuser . 


Recommendations 

The turbine stage is to be manufactured and tested. Hence, recom- 
mendations for further modification to improve stage performance will depend 
largely on the measured performance characteristics of the stage. Neverthe- 
less, a critical review of the design standard and predictions based on this 
standard lead to the following recommendations. 
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1. Two or more additional nozzle vane geometries should be de- 
signed and tested. This recommendation is based on the uncer- 
tainty of current prediction techniques for nozzle rows with 
axial passage width contraction. Since the nozzle row can 

be characterized by a number of geometric variables and aero- 
dynamic parameters, changes in the row solidity by a reduc- 
tion in the number of blades and by a reduction in blade 
chords (by reducing the i n 1 et-to-exi t diameter ratio) are 
recommended. 

2. Based on a review of the surface velocity distributions for 
the tandem row configuration with 16 blades in each row and 
boundary- 1 ayer calculations for the computed distributions, 
it would appear likely that further improvements could be 
made in the design-point performance of the rotor. it Is 
recommended that a design with more blades in the first row 
and fewer blades in the exducer section be considered. A 
design with 18 or 20 blades in the first row and 9 or 10 in 
the second would be recommended, subject to more detailed 
flow analyses of design alternatives. 
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TA8LE I 

MAJOR DESIGN GEOMETRIC VALUES FOR DATUM TURBINES 
USED IN THE PARAMETRIC INVESTIGATION 


Parameter 

Sub j ect 
Turbine 

Datum 3 

Datum 2 

Datum 1 

Reference 
Turbi ne 

Nozzle Vane Inlet 
Diameter, in (cm) 

7.654 
(19* 44 1 ) 

7.430 

(18.872) 

7. 144 
(18.146) 

6.869 

(17.447) 

6. 656 
(16.906) 

Nozzle Vane Ex! t 
Diameter, in (cm) 

5.426 

(13-782) 

5.438 

(13.812) 

5.438 

(13.812) 

5.331 

(13.541) 

5.118 

(13.000) 

Rotor 1 nlet Diame- 
ter, in (cm) 

5.268 

(13-381) 

5.281 

(13.414) 

5.281 

(13.414) 

5.177 

(13.149) 

4.970 

(12.624) 

Rotor Exi t Ti p 
Diameter, in (cm) 

3.460 

( 8 . 788 ) 

3.460 

(8.788) 

3.460 

( 8 . 788 ) 

3.460 

( 8 . 788 ) 

3.460 

( 8 . 788 ) 

Rotor Exit Hub 
Diameter, in (cm) 

1.822 

(4.628) 

1.822 

(4.628) 

1.822 

(4.628) 

1.822 

(4.628) 

1.822 

(4.628) 

Nozz 1 e Inlet Width, 
in (cm) 

0.560 

(1.422) 

0.434 
(1 . 102) 

0.417 

(1.059) 

0.401 

(1.018) 

0.385 

(0.978) 

Nozz 1 e Exit Wi dth , 
i n (cm) 

0.400 

(1.016) 

0.400 

(1.016) 

0.400 

(1.016) 

0.385 

(0.978) 

0.385 

(0.978) 

Rotor Inlet Wi dth , 
in (cm) 

0.400 

(1.016) 

0.400 

(1.016) 

0.400 

(1.016) 

0.385 

(0.978) 

0.385 

(0.978) 

Axial Length, in 
(cm) 

1.552 

(3.942) 

1.627 
(4. 132) 

1.694 

(4.303) 

1.765 

(4.483) 

1.765 

(4.483) 

Number of Rotor 
B 1 ades 

10 

1 1 

1 1 

1 1 

1 1 

Number of Sp 1 i tters 

10 

1 1 

1 1 

1 1 

1 1 

Number of Nozzle 

15 

15 

14 

13 

13 


Vanes 


The ''subject turbine 11 is essentially that selected at the conclusion 
of the preliminary design investigation. The rotor configuration was 
later changed to a tandem row. 

The "reference turbine" is the AiResearch design of Reference 2. 
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TABLE I 1 

SUMMARY OF STAGE PERFORMANCE AND 
OPERATING CONDITIONS 


Subject Turbi ne Reference Turbi ne 

(Rotor Diameter - (Rotor Diameter ~ 

Parameter 5.268 in) (13.381 cm) 4, 97 ? n) ( 12, 624 cm) 


Total -to-Total Efficiency 
at Rotor Exit, Ut* 


Total -to-Static Efficiency 
at Rotor Exit, T? 

'■'64- 

Total -to-Tota 1 Efficiency 
at Diffuser Exit/V^ 


Tota 1 -to-Stat i c Efficiency 
at Diffuser Exi t , 

SB * 

Over-All Efficiency, 


Inlet Total Temperature, 

Tqo ’ deg R 


I nlet Total Pressure, 
ps i a 

Weight Flow, TYl , 1 bm per 
sec 



Shaft Work, 



k i 1 owa 1 1 s 


Turbine Speed, fs/ , rpm 

. . 3/4 

Specific Speed, f\L , ft / 
(mi n) (sec 1/2) S 


0.9433 

0.9009 

0.9358 

0.9300 

0.9249 

2060 (1144 deg K) 

25 (172 kN/m 2 ) 

0.7484 (0.3395 kg/sec) 

18.08 

36,000 

3/4 

75-1 (30.8 rn / 

(mi n) (sec 1 /2) 


Turbine Inlet Total-to- 
Diffuser Exit Static Pres- 
sure Ratio, D 

rn 


00/ p P 


1.763 


0.913 

0.870 

0.896 

0.890 

0.887 

2060 (11 44 deg K) 

25 (172 kN/m 2 ) 

0.7484 (0.3395 kg/sec) 

17.33 

36,000 

75.9 (31.1 m 3/4 / 

(mi n) (sec 1/2) 

1.763 




Over-all efficiency ( f ( n s/) ~ ('total -to-static efficiency at diffuser 
exit; - efficiency decrement due to disk friction). 



REPRODUCIBILITY OP THE: ORIGINAL PAGE IS POOR 
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TABLE I I 1 

PRELIMINARY DES I GN**P0 1 NT PERFORMANCE AND FLOW FIELD DETAILS 

UcbxwtS uPLuATiNo CoNJlTiONS 


ROTATIONAL SPEED HP A ’“‘36009. GO 


..AoSOLOlfc. MASS FLO „ W _ • 7 4 b V 

"TOTal-TO-TOTAU PHfcSSUKt .W^Tin a‘T*V PkT 4 " X * V *- V 5i 

To i |i» L“ 1 1 7 a C P P L G 3 ONE K' A T 1 0 a T 4 . PA'S 4 i * 7 v 9 7 

ToTAL-To-TOTAL>KLSSUhJi RATIO AT S PhTg *"l,7»60 

.TOTAL-1 0 " *> I /■. T 1 C . ■ p R £ S S U R £ , W A T 1 0 . A L > P R S t* 1 * 7 1> 3 0 

G A S C 0 J\ S I a N T R G A S A b • 4 j 7 0 

.SPECIFIC ntAl HaTIO QfiA J. 1,6667. 


TURalMt PlYnOluS NUhgjtK RfcYT 7 2 u 4 1 

.SP£C 1 f 1 C . SP & fc U. SP So.. 76*10 

m 


STaGE PERFORMANCE 


jSH AF T .WpRK . 

’overall LfFYcltNCY '* 

JOTal-TO-T OT a l - fcF f I c 1 1 N C V a T 4 ... 
T0 T A L-TQ-ST a T 1C cFKj.clfc.NCY A T A 
TOTAL-TO-TOTAU EFFlClENCy AT 6 
T0 Tal-T0-STaTaC efficiency at s 


work 

16,077. 

etaov 

• 9249 

ETAT 4 

.*9433 

etasa 

• 9 GG 9 

ET A Tb 

,935b 

ETASb 

•9300 


COMPONENT EFFICIENCY DECREMENTS 


..NOZZLE TpT At...; detan '. ,0161 

PQTOH TOTAL , DETaP *04U2 

. OlFF^SfcR J0T A L OE t AET 3 7 „ 


6 K fcA K U 0 « N *pF “E p p i c I E N C y’D E C R E M E N T $" 


scroll f RxcT Ton ~ ” ‘detasc “ ~ .0026 

_ NOZZLE. VAnE LOAL/JnG _ L OETaVb « 0062 

NOZZLE aSPLCI R a V 1 0 OlTaVa * 0033 ' 

NOZZLE VANE TRAILING EDGE DE.TaVT ^_.-.0002 

NOZZlz K E Y N g L u S NUHbER 1 DE^ AVP “*GuGv 

. V aNcLES S SP;.cE PRIvllON • deTaV’L ,0032 

POTOP r R iCT I ON " "DcTaRF ,o079 

ROTOR CLEARANCE. . DETaKC ,*0109. 

PoTOH DIFFUSION AND FRICTION OeTa^D •0£9<t 

. DIFFOStP 1 GE CWl 1,0077 

DlSCrtAKGE UETaOE .0060 

DISC FkICTIoN CWINDAGt) DETaOF *00bl 

NOZZLE VaNE ToTal ^ DETaVN “',0103* 


TABLE 111 (CONTINUED) 

PRELIMINARY DESIGN-POINT PERFORMANCE AND FLOW FIELD DETAILS 
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FLO A MO GEOMETRIC LA»* mT STATION 0 


AsSOLUlE TOTAL FISSURE POO 3600.00 

AaSOLUTg ’luiAu'TEMHfcHATURe"*”” TqO " 2060*00 

STATIC PRESSURE . PO ' 3595*34 

STATIC Tt -ttHtHATupe TO ?0btt«93 

AbSOuolc VLLOClTy Co 3 6,23 

ABSOLUTE i v 'AC« NUMBER EMQ *0394 

M£aN ^CmOLL) . . dOmjn l6*o*0 

MEAN MyUKAULiC UIAMETER ( Sc«OLL) DHSCIN 2,546 

scroll friction >. coefficient _ cf.sc *oo&3i 


FLOW ANO GfiOrtfiTWIc. DaTA At STATION 1 


AosOLUTt TOT^l PR’ESSUrS 

absolute total temperature 
static pressure 

STATIC TENPeKaTuHE 
■AB SOLUTE VtLOciTY 
M^RlijjuN^L velocity 

ABSOLUTE tangential velocity 

abSolu i E GAS A^OLE . ... 

absolute mach number 
OIAM&TER 

PaSSaGE w 1 0 T rt 


POl 

3595 * 00 

* 0 i 

2060*00 

Pi 

356 1*92 

Ti 

20b2. 0 

Cl 

ibO.ii 

CM 1 

_ 65*02 

CT1 

123,71 

ALFjIN 

- 55-50 

EM 1 

• 1 0 £> 

DliN_... 

.7,65a 

BllN 

• 56Q0 





REPRODUCIBILITY OF THE ORIGINAL PAGE. IS POOR. 
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TABLE II! (CONTINUED) 

PRELIMINARY DESIGN-POINT PERFORMANCE AND FLOW FIELD DETAILS 
A\U G E 0 M £ T fi I C, . 0 A T A A T S T A T I On 2 
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ABSOLUTE TOTAL PRESSURE 

absolute total temperature 

STATIC PRESSURE 
STATIC TEmPERaTURE 
ABSOLUTE VELOCITY 

meridional' velocity 
ABSOLUTE T mnGENY i a l velocity" 
ABSOLUTE Gas ANGLE _____ 

absolute Mach mumper 

diameter . . __ 

Passage width ’ * 

JVO Z 2 L £_ JK.A I L 1 IMG £ Dpt.. A NGLE . _ 

NOZZLE TRAIL I NB EDGE THICKNESS 
-NUMBER. .OF NOZZLE .VaN£S 
N.022LE ThRO A T/HtIGh1 ASPECT RaTIO 

..R'aXI 9~Q£.JUSa I L ji NCj..E06E./J.HR.0aT... ^ I pT H.. 


P02 

3b7b,19 

Tf;2 

.. 2 0 6 U • 0 0 

P2‘ 

2 t>$ 6*63 

T2 

_ .1693*67 

C2 

’702.0a 

CM2 

184*96 

CT2 

677-27 

ALP’S IN 

7 4 * 7-3 

£M2 . 

* 5 1 3 

02 IN 

5 . 426 

b2In 

. 4 j 2 0 

BETAS.!... 

74. 2.3 

T£Nin 

’0100 

Zjvj 

.is 

AR 

1 o623 

_ XeOTW. 

0 033 


flow and geometric daTa "at" 

station 3" 



.ABSOLUTE JOTaL pressure 
absolute total temperature" 

_R£LATIVe.._TOT a L PRESSURE 

P03 

To3 

P03R 

: 3569,66 

2060c 00 
• 2925.42 

Relative tot A l Tehper a tur £ 

^STATIC PRESSURE 

T 03R 
P3 

■ iy02-31 
2652 64 

STATIC TEMPERATURE 

absolute velocity 

T3 

_ _C3 

3.383 ® 23 
‘723-76 

meridional velocity 

absolu i e.JaN6enTi a l velocity 

CM3 

CT3 

1 9 8 ’ 1 5 
696°11 

absolute g A s a hgle 
absolute hach number 

ALF3IN ' 
EM3 

74-11 
°53 1 

relative velocity 
Relative Tangential velocity 

W3 

_wt3 

237.75’ 

-l3l.3q 

relative gas angle’ 
Relative mach number 
'WHEELS PE Eu 
DIAMETfcR 

“aLF3«Y~ 
em 3 r 
U3" " 

D3IN 

"-33.5b 
-174 
8*2 7 - 49 
5-268 

passage width ✓ 

_RO_l OR inleT blade ANGLE 

83 1 n 

JVE.TA3I 

*4 00 0 
4,78 

rotor blade number 

ZR3 

12 


SI 


LV V. - J X •? L xii **. 


-..s* £ t i.tti,VrJi 




TABLE 1H (CONTINUED) 

PRELIMINARY DESIGN-POINT PERFORMANCE AND FLOW FIELD DETAILS 


FLOW AND. GEOfCTKIC OaTa aT STATION 4 


total absolute . pressure po4 <soe>o # 

TOTAL AbSuLUTE TEMPERATURE " , Tq4 1671*24 

TOTAL KP.LaT 1 VE PRESSURE _J Pfl4R . . 2244*10 

Total relative TekPeraTure to4r 1729.31 

STATIC PRESSURE _ __ . R4 200‘0.33. 

STATIC TEMPeRa JUKE T4 ioibl.b7 

ABSOLUTE" VELOCITY C4 _ 241*42 

MERIDIONAL VELOCITY ~ CMrt4. 241*42 

absolute ta-oe^Ti a l velocity ctm4 • o„gg 

ABSOLUTE GAS a^oLe ' ALFH4I ‘ 0.00 

absolute mach number em4 ,109 

Relative velocity ' w 4 " 479*93 

RcLATivt -Ia^entiaL velocity _ *tm 4 _..~4i4*6s 

Relative gas a'^gle ALF 4 RI "59*60 

RE LA r I vE MACH (NyMBhR £M4R „3?6 

Wh£ELSPEED w UM 4 4 14 9 85 

DIAMETER aT^HUB _ DH4.IN, 1*022. 

diameter aJ m eaN Radius DM 41 N 2°o4i . 

Dl.A«ETtK A t . ShRUUp; DS4IN 3*460 

PaSSaGe width ~ ” B4lN ' ” 00 I 9 O 

•_ 8 Lade ea.it a ngle a t mean r a diu$ „.6&‘Ta4I_ " 56 * 07 * 

TRAILING edge ThIcRHeSS' aT mean 'RadIUS“'Te«I^' ~ <. 0 b 00 ; 

MEAN KOTOR BLADE PASSAGE JTYDRAULlC_DlA_ DHRI.nl _ g^GoSQ 

ROTOR'fiLAOE NUMBER '• * 2R4 ~ **12 





REPRODUCIBILITY OF THE ORtGINAL PAGE IS POOR. 
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TABLE Hi (CONTINUED) 

PRELIMINARY DESIGN-POINT PERFORMANCE AND FLOW FIELD DETAILS 


' K LW" A ND^eeOMET « I c DaTa at station s 
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ABSOLUTE T*OtAL PRESSURE * 

absolute total temperature 
static pressure 

STATIC T E M P E H A T 0 K’ [£ 

ABSOLUTE VELOCITY 

ABSOLUTE hAC.ri. NUMBER 

ShROUq d 
HUB OiAMElEw 

Dl’FF US l« a^EA *A t IO-OV'eR a l‘l 

U i F r US^K A H L A ^aT - lO* , ' , EQU < ^Vq l 
DIFFUSER ARfcA RATIO-EqUiVo i 

.D XFFUSeR,„^UrHaLI LE^G t Fi«*( 



I X r r UULH I I^LULJVHwU 

OIFFuSER lO$S COEFFICIENT 
DIFFUSER aVeRaBE-TO-MeaN LOSS Ra’TIU 


pos 

■<£050. 11 

Tob 

it 71*24 

p 5 

2 0 o, I <.97 

T S 

1 6 o 6 » Sb 

C5 

8 b *70 

EMb 

o 069 

DSbIN 

£> ° 1 60 

_ .DH^IN 

1 & 

aRdif 

2 » 69 

ARDIFU 

2,69 

APOIFL 

2*69 

JOTh_ELunL. 

8*75. 

:a eldna 

11*26 

ELDI'NL 

7 « 1 6 5 

ELUINA 

9 0 b4 ^ 

alf'dil 

.5 *S3 

ALFoIa 

4*29 

DTFP'F 

0 799 

8LK4 

0 1 0 0 

CLOSDA 

Q 173 

CLOSqR 

1 p 0 0 0 


55 


TABLE IV 

MERIDIONAL FLOW PATH SPECIFICATION 
OF THE BRU ROTOR 


Fi rst Blade Row 

Outer Contour Coordinates'' Hub Contour Coordinates 


Radi us 
(in) (cm) 

2.634 6.690 

2.510 6.375 

2.397 6.088 

2.293 5.824 

2.197 5.580 

2.105 5.347 

2.018 5.126 

1.938 4.922 

1.835 4.661 


1.835 4.661 
1.750 4.445 
1.733 4.402 
1.730 4.394 
1.730 4.394 


Axial Distance 
Measured from 
Rotor Exi t 

Hill i£2 1 

- 1.152 - 2.926 

-1.141 - 2.898 

-1.119 -2.842 

-1.085 -2.756 

- 1 . 04 1 -2 . 644 

-O .983 -2.497 

-0.913 -2.319 

- 0.830 - 2.108 

- 0.678 - 1.722 


Radi us 
(in) (cm) 

2.634 6.690 
2.427 6.165 
2.230 5.664 
2.039 5.179 
1.853 4.707 
1.669 4.239 
1.493 3.792 
1.329 3.376 
1.122 2.850 


Axial Distance 

Measured 

from 

Rotor 

Exit 

Iij2l 

(cm) 

-1.552 

-3 . 942 

-1.516 

-3.851 

-1.469 

-3.731 

-1.411 

-3.584 

-1.338 

-3.398 

- 1 . 244 

-3.160 

-1 . 120 

-2.845 

-O. 96 I 

-2.441 

-0.678 

-1 .722 



Second 

Blade Row 




-0.678 

-I .722 

1 . 1 22 

2.850 

- 0.678 

-1.722 

-0.451 

-1.145 

1.018 

2.586 

-0.475 

- 1.206 

-0.313 

-0.795 

0.955 

2.426 

-0.297 

-0.754 

-0.198 

-0.503 

0.925 

2.349 

-0. 148 

-O .376 

-0.059 

- 0.150 

0.911 

2.314 

+0.018 

+0.046 


Rotor casing coordinates are obtained with the addition of a 0.010 in 
(0.0254 cm) c 1 ea ranee. 
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TABLE V 

FIRST BLADE ROW DEFINITION ALONG THE STREAMTUBE CENTERS 


Outer Casing Streamtube 
Axial Distance 




Measured 

1 from 

Angular 

Normal 

Radius 

Rotor 

Exit 

Coordinate 

Thickness 

(in) 

(cm) 

.Onf^ 

(cm) 

,( de 9). 

(in) 

IcmL 

2.634 

6.690 

-1.170 

-2.972 

0.00 

0.0500 

0. :l 270 

2.568 

6.523 

-1.163 

-2.954 

- 0.37 

0.0500 

0.1270 

2.506 

6.365 

-1.155 

-2.934 

- 0.95 

0.0500 

0.1270 

2.446 

6.213 

-1.146 

-2.911 

- 1.65 

0.0500 

0. 1270 

2.389 

6.068 

-1.134 

-2.880 

- 2.45 

0.0500 

0. 1270 

2.333 

5.926 

-1.119 

-2.842 

- 3.33 

0.0500 

0. 1270 

2.279 

5.789 

-1.102 

-2.799 

- 4.25 

0.0500 

0.1270 

2.227 

5.657 

-1.082 

-2.748 

- 5.23 

0.0500 

0.1270 

2. 178 

5.532 

-1.056 

-2.682 

- 6.24 

0.0500 

0. 1270 

2.131 

5-413 

-1.026 

-2.606 

- 7.27 

0.0500 

0. 1270 

2.086 

5.298 

-0.993 

-2.522 

- 8.31 

0.0500 

0. 1270 

2.042 

5.187 

-0.958 

-2.433 

- 9.36 

0. 0500 

0. 1270 

1.999 

5-077 

-0.920 

-2.337 

-10.39 

0.0500 

0. 1270 

1.958 

4.973 

-0.879 

-2.233 

-11.38 

0.0499 

0. 1267 

1.918 

4.872 

-0.834 

-2.118 

-12.33 

0.0482 

0. 1224 

1.880 

4.775 

-O.786 

-1.996 

-13.21 

0.0436 

0.1107 

1.846 

4.689 

-0.733 

-1.862 

-14.00 

0.0359 

0.0912 

1.815 

4.610 

-0.678 

-1.722 

-14.69 

0.0250 

0.0635 




Hub Streamtube 



2.634 

6.690 

-1.538 

-3.906 

1.68 

0.1000 

0.2540 

2.530 

6.426 

-1.518 

-3.856 

.1.06 

0.1016 

0.2581 

2.430 

6. 172 

-1.497 

-3.802 

0.13 

0 . 1 03 1 

0.2619 

2.334 

5.928 

-1.474 

-3. 744 

- 0.98 

0.1046 

0.2657 

2.239 

5.687 

-1.449 

-3.680 

- 2.20 

0.1061 

0.2695 

2. 145 

5.448 

-1.422 

-3.612 

- 3-51 

0. 1076 

0.2733 

2.053 

5.215 

-1.393 

-3.538 

- 4.89 

0.1091 

0.2771 

1.961 

4.981 

-1.359 

-3.452 

- 6.31 

0. 1106 

0.2809 

1.871 

4.752 

-1.321 

-3.355 

- 7.77 

0.1121 

0.2847 

1.782 

4.526 

-1.278 

-3.246 

- 9.22 

0.1136 

0.2885 

1.694 

4.303 

-1.228 

-3.119 

-10.63 

0. 1152 

0.2926 

1 .610 

4. 089 

-1 . 171 

-2. 974 

-11.96 

0.11 67 

0.2964 

1.529 

3.884 

-1.105 

-2.807 

-13.16 

0.1183 

0.3004 

1.454 

3.693 

-1.030 

-2.616 

-14.15 

0.1200 

0.3048 

1.382 

3.5IO 

-0.949 

-2.410 

-14.86 

0.1146 

0.2911 

1.316 

3.343 

-0 .862 

-2.189 

-15.21 

0.0970 

0.2464 

1.255 

3. 188 

-0.772 

-I.96I 

-15.16 

0.0671 

0.1704 

1.200 

3.048 

-0.678 

-I.722 

-14.70 

0.0250 

O.0635 
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TABLE V) 

SECOND BLADE ROW DEFINITION ALONG THE STREAMTUBE CENTERS 


Outer Casing Streamtube 


Radi us 

Axial Distance 
Measured from 
Rotor Exit 

Angular 
Coord i nate" 

Normal 

Thickness 

On) 

(cm) 

(in) 

(cm) 

(deg) 

.(in). 

(cm) 

1.815 

4.610 

-0,678 

-1.722 

10.255 

0.020 

0.051 

1.785 

4.534 

-0.620 

-1.575 

9.337 

0.061 

0.155 

I.76I 

4.473 

-0.561 

-1.425 

8.114 

0.071 

0. 180 

1.741 

4.422 

-0.503 

-1.278 

6.628 

0.070 

0.178 

1.726 

4.384 

-0.446 

-1.133 

4,806 

0.068 

0.173 

1-715 

4.356 

-0.393 

-0.998 

2.709 

0.061 

0.155 

1.708 

4.338 

-0.341 

-0.866 

■0.399 

O.O56 

0. 142 

1 . 705 

4.331 

-0.293 

-0.744 

- 2.096 

0.048 

0. 122 

1.703 

4.326 

-0.246 

-o. 625 

- 4.710 

0.044 

0.112 

1.701 

4.321 

-0.201 

-0.511 

- 7.400 

0.043 

0.109 

1.700 

4.318 

-0.157 

-0.399 

-10.135 

0.038 

0.097 

1.700 

4.318 

-0. 1 1 5 

-0.292 

-12.917 

0.034 

0.086 

I.700 

4.318 

-0.059 

~o. 150 

-16.700 

0.026 

0. 066 

1.179 

2.995 

-0. 678 

Hub Streamtube 
-1.722 2.342* 

0.025 

0.064 

1.150 

2.921 

-0.622 

-1.580 

2.867 

0.083 

0.21 1 

1.124 

2.855 

-0.565 

-1.435 

3.315 

0. 122 

0.310 

1 .099 

2.791 

-0.507 

-1.288 

3.389 

0. 147 

0.373 

1.076 

2.733 

-0.448 

-1.138 

2.989 

0.154 

0.391 

I.056 

2.682 

-0.389 

-0.988 

2.219 

0. 147 

0.373 

1.037 

2.634 

-O.330 

-0.838 

1.136 

0.131 

0.333 

1.020 

2.591 

-O.272 

-0.691 

- 0.338 

0.115 

0.292 

1 . 006 

2.555 

-0.216 

-0.549 

- 2.384 

0.096 

0.244 

0.996 

2.530 

-0.161 

-0.409 

- 4.833 

0.076 

0.193 

0.989 

2.512 

-0. 108 

-0.274 

- 7. 565 

0.056 

0. 142 

0.984 

2.499 

-0.057 

-0,145 

-10.554 

0.039 

0.099 

0.978 

2.484 

-0.008 

-0.020 

-14.849 

0.025 

0.064 


-I- ' . ' : 

“Angular coordinates are measured with reference to the stacking point 
of the second blade row axial blade sections. 
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TABLE VII 


Radius 

NOZZLE 

VANE GEOMETRY 

Angular 
Coordi nate 

SPECIFICATION 

Tangential 

Thickness 

(in) 

(cm) 

.(, de .g). 

(in) 


3.827 

9.720 

0.00 

0.0102 

0.0259 

3.774 

9.586 

1.09 

0.1300 

0.3302 

3.720 

9.449 

2.20 

0.1769 

0.4493 

3.665 

9.309 

3.31 

0.2082 

0.5288 

3.610 

9.169 

4.43 

0.2312 

0.5872 

3 . 555 

9.030 

5.56 

0.2482 

0.6304 

3.499 

8.887 

6.72 

0.2615 

0.6642 

3 « 445 

8.750 

7.91 

0.2715 

0.6896 

3.392 

8.616 

9.13 

0.2786 

0.7076 

3.341 

8.486 

10.39 

0.2834 

0.7198 

3.292 

8.362 

11.70 

0.2871 

0.7292 

3.246 

8.245 

13.06 

0.2894 

0.7351 

3.204 

8. 138 

14.47 

0.2896 

0.7356 

3.163 

8.034 

15.92 

0.2875 

0.7302 

3.126 

7.940 

17-42 

0.2822 

0.7168 

3.091 

7-851 

18.95 

0.2735 

0.6947 

3.057 

7.765 

20.51 

0.2611 

0.6632 

3.026 

7.686 

22.10 

0.2452 

0.6228 

2.996 

7.610 

23.72 

0.2263 

0.5748 

2.967 

7.536 

25.36 

0.2045 

0.5194 

2.939 

7.465 

27.02 

0.1818 

0.4618 

2.912 

7.396 

28.70 

0.1603 

0.4072 

2.885 

7.328 

30.40 

0.1405 

0.3569 

2.859 

7.262 

32.13 

0.1224 

0.3109 

2.834 

7.198 

33.87 

0.1063 

0.2700 

2.809 

7.135 

35.63 

0.0924 

0.2347 

2.784 

7.071 

37-40 

0.0806 

0.2047 

2.760 

7.010 

39-20 

0.0710 

0. 1803 

2.737 

6.952 

41.01 

0.0634 

0. 1610 

2.713 

6.891 

42.84 

0.0577 

0. 1465 
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TABLE VIII 

DIFFUSER SPECIFICATION 


Axial 

Di stance 

from 

Rotor Exit 

,(tn), 

(cm) 

0.00 

0.00 

1.00 

2.54 

1.25 

3.17 

2.00 

5.08 

3.00 

7.62 

4.00 

10.16 

5.00 

12.70 

6.00 

15.24 

7.00 

17.78 

8.25 

20.95 

9.25 

23.49 


Casi ncj 

Radius 

,( i n) 

(cm) 

1.730 

4.394 

1.736 

4.409 

1.740 

4.420 

1 . 786 

4.486 

1.815 

4,610 

1.875 

4. 762 

1.950 

4.953 

2.047 

5.200 

2.202 

5-593 

2.500 

6.350 

2.580 

6.553 


Hub Radius is Constant 
for the Entire Length = 
0.91 1 in (2.314 cm) 


Note : Area schedule is based on an assumed constant static pressure be- 

tween 0 and 1.0 in (2,54 cm), linear static pressure rise between 
1.0 in (2.54 cm) and 8.25 in (20. 95 cm), and a linear drop of total 
pressure with axial distance from 0 to 8.25 in (20.95 cm). 



TABLE IX 


TURBINE INLET SCROLL GEOMETRY ‘SPEC IF I CAT ION 


Maximum Radius 

Azimuth Angle 


Radius of 

of the Cross 

Measured from 

Area of the 

the Cross 


Section, Rm 

Scroll Inlet 

Cross 

Section 

Section, 

R0R 


R0L 


list 

(cm) 

(deg) 

(sq In) 

(sq cm) 

On) 

(cm) 

(In) 

(cm) 

(In) 

(cm) 

8,427 

21.404 

0.0 

14.209 

91.671 

2.050 

5.207 





8.301 

21.084 

20.0 

13.372 

86.271 

1-987 

5.047 





8.171 

20.754 

40.0 

12.537 

80.884 

1.922 

4.882 





8.038 

20.416 

60,0 

11.704 

75.510 

1.855 

4.712 





7.899 

20.063 

80.0 

10.875 

70.161 

1.786 

4.536 





7.756 

19.700 

100.0 

10.049 

64.832 

1.715 

4.356 





7.608 

19.324 

120.0 

9.226 

59 . 522 

1.641 

4.168 





7-453 

18.930 

140.0 

8.408 

54.245 

1.563 

3.970 





7.291 

18.519 

160.0 

7.594 

48.993 

1.482 

3.764 





7.120 

18.085 

180.0 

6.785 

43.774 

1.397 

3.548 





6.940 

17-628 

200.0 

5.983 

38.600 

1.307 

3.320 





6.747 

17.137 

220. 0 

5.186 

33.458 

1.211 

3.076 





6.536 

16.601 

240.0 

4.395 

28.355 

1.108 

2.814 





6.301 

16.004 

260.0 

3.611 

23.297 

0.998 

2.535 





6.033 

15.324 

280.0 

2.838 

18.310 

0.876 

2.225 





5«718 

14.524 

300.0 

2.079 

13.413 

0.741 

1.882 





5.326 

13.528 

320.0 

1.341 

8 . 65 1 

0.587 

1.491 





4.784 

12.151 

340.0 

0.638 

4.116 

0.408 

1.036 





4.394 

11.161 

350.0 

0.306 

1-974 

0.315 

0.800 





4.301 

10.924 

352.0 

0.242 

1.561 

0.300 

0.762 





4.203 

10. 676 

354.0 

0. 180 

1.161 

0.287 

0.729 





4.152 

10.546 

355.0 

0.149 

0.961 

0.282 

0.716 





4.084 

10.373 

356.0 

0.119 

0.768 



0.227 

0.576 

0.258 

0.655 

4.006 

10.175 

357.0 

O.O 89 

0.574 



0. 1 58 

0.401 

0.179 

0.455 

3.940 

10.008 

358.0 

0.059 

0.381 



0.100 

0.254 

0.113 

0.287 

3.881 

9.858 

359.0 

0..029 

0.187 



0.048 

0. 122 

0.054 

0.137 

3.827 

9.720 

360.0 

0.000 

0.000 



0.000 

0.000 

0.000 

0.000 

Note; 

The nomenclature of the cross 

sections 

is defined 

in Figure 

22. 




a\ 
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TABLE X 

ROTOR ADIABATIC WALL TEMPERATURE AND BIOT NUMBER 
DISTRIBUTION ASSUMED FOR ROTOR TEMPERATURE 
AND STRESS CALCULATIONS 




Front Face of Disk 



Back Face of 

Disk 

.. * 
Rad t us 

Ratio 

•/ride 

_, :t Temperature 
Region' Ratio 

Biot 

Number 

Region 

Temperature 

Ratio 

Biot 

Number 

0,000 

EF 

0.8032 

0.000 

AB 

0.3906 

1 . 505 

0.076 


0.8030 

0.065 


0.3906 

1.505 

0.114 


0.8025 

0.084 

BC 

0.9631 

0.084 

0.190 


0.8013 

0.113 


0.9636 

0.113 

0.266 


0,7395 

0.139 


0.9642: 

0.139 

0.304 


0.7983 

0.150 


0.9646 

0.150 

0.342 


.0.7971 

0.161 


0.9651 

0..161 

0.346 

DE 

0.8203 

3.500 


0.9656 

0.162 

0,418 


0.8369 

2.040 


0.9667 

0. 182 

0.494 


0.8476 

1.320 


0.9674 

0.201 

0.570 


0.8579 

1.180 


0.9689 

0.219 

0. 646 


0.8682 

1,050 


0.9706 

0.237 

0. 661 


0.8706 

1,040 


0.9714 

0.240 

0.677 


0.8725 

1.040 


0.9714 

0.243 

0.692 


0.8735 

1 . 050 

BC 

0.9716 

0,245 

0.722 


0,8784 

1.070 


0.9726 

0.253 

0.79S 


0.8901 

1 . 1 10 


0.9747 

0.268 

0.874 


0.9028 

1 . 060 


0.9772 

0.283 

0.950 


0.916 

1 , 200 

s 

0.9797 

0.298 

1.000 
J - 

CD 

1.000 

0.307 

CD 

1.0000 

0.307 


Radii are ratioed to the rotor inlet radius of 2.634 in -(6.690 cm). 
For definitions of regions, see Figure 35- 

Temperatures are normalized to the rotor inlet disk face temperature 
of 2048 deg R (1140 deg l<) , 



TABLE XI 


PHYSICAL PROPERTIES 


HIGH TEMPERATUF 


Test 
Tem- 
pera- 
ture 
(deg F) 


Stress 

Rupture 

(50,000 

hrs) 


0.4x10-5 

Per 

Cent/hr 
Minimum 
Creep Rate 


Endur- 
ance 
LI mi t 

TT58 

Cyc 1 es) 


Ultimate 

Strength 


Yield 
2 Per 
Cent 
Offset 


Casting 

Alloys 


7I3C 1400 40,000 27,000 26,000 135,000 109,000 -0.05 

(1033 
deg K) 


713LC 1400 

(1033 
deg K) 


INI 00 1400 

(1033 

deg K) 


43,000 22,000 


50,000 


25,500 138,000 


11,000 155,000 


110,000 - 0.10 


125,000 - 0.80 


35,000 



6 10 0.65 Developed to 

combine good 
stress-rupture 
properties wi th 
excellent room 
temperature, 
ductility, and 
strength; for 
use with Inte- 
grally cast 
turbine wheels 

10 19 0.59 Developed to op- 

timize stress 
rupture and room 
temperature, 
ducti 1 i ty , and 
strength 

6 7-2 0.22 Good strength 

to density char- 
acter? sties, 
long-time em- 
brittlement 
problem 


ON 

ro 



TABLE Xj (CONTINUED) 


PHYSICAL PROPERTIES OF HIGH TEMPERATURE ALLOYS 


A1 Toy 
Name 

Test 
Tem- 
pera- 
ture 
(deg F) 

Stress 

Rupture 

( 50,000 

hrs) 

Creep 
a. 4x10-5 
Per 

Cent/hr 
Mi ni mum 
Creep Rate 

Endur- 
ance 
Limi t 
(10 8- 
Cycles) 

U1 timate 
Strength 

UDIMET 

500 

1400 

(1033 

deg K) 

33,000 

25,000 

10,000 

120,000 

MAR- 

::200 

1400 

(1033 

deg K) 

30,000 

- . 

25,000 

135,000 


Wrought 

Alloys 


UDIMET 

700 

l4oo 
(1033 
deg K) 

43,700 


26,000 

145,000 

WASPAL0Y 

1400 

(1033 

deg K) 

24,300 


18,000* 

120,000 

RENE 41 

1400 
(1033 
deg K) 

26,800 

- 

15,000* 

100,000 


Assumes 6 -0.15* 

U( 


Yield 
2 Per 
Cent 
Offset 

Per 
Cent 
Wei ght 
Change 
(200 hrs) 

Per 

Cent 

Elonga- 

tion 

Per 
Cent 
Reduc- 
tion 
in Area 

Endur- 

ance 

Limit 

Stress/ 

Rupture 

Strength 

Remarks 

105,000 

-0.15 

21 

24 

0.30 

Good oxidation 
resistance, rela 
tively low 
strength 

122,000 

1 

o 

* 

Ml 

o 

3.5 


0.83 

Lacks ductility 
at. lower tem- 
peratures as 
needed for in- 
tegrally cast 
turbine wheels 

120,000 


8 

10 

- 

Difficult and 
expens i ve to 
forge 

100,000 


22 

30 


Properties not 
as good as 
casting alloy 
at 1400 deg F 
(1033 deg K) 

82,000 

- 

7 

15' 


Properties not 


as good as cast- 
ing alloy at 
1400 deg F 
(1033 deg K) 

o\ 

XjO 



TABLE XI (CONTINUED) 

PHYSICAL PROPERTIES OF HIGH TEMPERATURE ALLOYS 



Test 

Tern- 

pera- 

Stress 

Rupture 

Creep 

0.4xl0~5 

Per 

Cent/hr 

Endur- 

ance 

Limit 


Yield 
2 Per 

Per 
Cent 
Wei ght 

Per 

Cent 

Per 

Cent 

Reduc- 

Endur- 
ance . 
LImi t 
Stress/ 


Alloy 

ture 

(50,000 

Minimum 

(10*“ 

Ultimate 

Cent 

Change 

Elonga- 

tion 

Rupture 


Name 

(deg F) 

hrs) 

Creep Rate 

Cyc 1 es) 

Strength 

Offset 

(200 hrs) 

tion 

i n Area 

Strength 

Remarks 

INCONEL 

718 

1400 

(1033 

deg K) 

20,000 


12,000* 

80,000 

62,000 


10 

25 


Insufficient 
strength at 
1400 deg F 
(1033 deg K) 
difficult to 
machine 


Assumes 




/ 




« 0* 15> 


G\ 



TABLE XII 


SUMMARY OF DISK 

STRESSES AND 

DEFLECTIONS 



Steel (Cold) 


Steel (Hot) 


Aluminum (Cold) 

120 Per 

100 Per 

120 Per 

L r 

140 Per 

150 Per 


Property 

Cent Speed 

Cent Speed 

Cent Speed 

Cent Speed 

Cent Speed 

Average Tangential Stress (psi) 

12,717 

12,258 

18,448 

25,757 

4,861 

(kN/m 2 ) 

87,680 

84,5*6 

127,194 

177,588 

33,515 

Axial Deflection at Tip (in) 

+0.0008 

0.0016 

0.0019 

0.0023 

0.0005 

(cm) 

0.0020 

0.0041 

0.0048 

0.0058 

0.0013 

Radial Deflection at Tip (in) 

+0.0013 

0.0398 

0.0402 

0.0407 

0.0007 

(cm) 

0.0033 

0.1011 

0. 1021 

0. 1034 

0.0018 


c\ 

vn 



66 


TABLE XIII 


SUMMARY OF BLADE' FREQUENCY AND VIBRATORY STRESS RESULTS 


Blade 

Row 

Location 

Method 

(cps) 

X (E) 

K 

**"'i 

b 

(psi) 

_V 

(psi) 

1 

L.E. 

Flat Beam (1) 

13,450 

22.4 

1.3 

2.5 

110 

358 

1 

T.E. 

Voysey (2) 

31,800 

53.0 

1.3 

2.5 

0 

0 

2 

L.E. 

Cant i lever 
Plate (3) 

3,940 

6.8 

1.2 

6.0 

1 485 

10,620 

2 

T.E. 

Canti lever 
Plate (3) 

3,940 

6.8 

1.2 

6.0 

1260 

9,090 


(1) Mechanical Design and System Handbook - Rothbart - pp. 6-58 


( 2 ) 


(3) 


4- « C y. Y'" v \ (y K 

T n 11 jy. 


for Flat Beam.C^ = 11.30 


Some Vibration Problems in Gas Turbine Engines - Voysey 

for Trailing Edge Panel 


fi 


y\ 


h G A J 

J- f\ G 


21. 


Where 


i 4TTL2- 

L. - Diagonal Length a 


E U, 


.3 5 


=35 degrees from Reference Through 
Blade Tip Parallel to Center Line 

Mechanical Design and System Handbook - Rothbart - pp. 6~59 

4- , 

for Cantilever P late, C*<yi = 3 • 3 1 


-f = JL xlo X k 


on 


71 a 2 - , 

where f or T\ =1, p\ = Blade Thickness « 

(X = Blade Height = JL 
= B 1 ade Wi dth 

1^^= Material Constant - 1.0 for Steel 
Y" - Radius of Gyration = in 

£» = Young's Modulus, psi 

^ « Material Density, lb/in3 

= Constant (a Function of Taper Ratio) 


The vibratory stress is calculated from the equat ion^v- KfcK where 

= bending stress, K = stress concentration factor, and 0< = ampli- 
fication factor (Ref 15) given in the above table. 
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TABLE XIV 

TURBINE RESEARCH PACKAGE ASSEMBLY AND TOOL LIST 


Part Number Description Quantity 


1 1 59-D004 



Wheel Turbine 

1 

1 159-E005 



Wheel Assembly Turbine 

1 

1159-R006 



Scroll Fabrication Assembly 

1 

1 1 59-RQ07 



Scroll Machining Assembly 

1 

1 159-R008 



Nozzle Machining Assembly 

1 

1 1 59-D009 



Shroud Machining 

1 

1 159-DO 10 



Diffuser Machining 

1 

1 1 59-DO 1 1 



Exhaust Adapter 

1 

1159-0012 



Rotor Back Shield 

1 

11 59-BO 13 



Pin Cap 

5 

1 1 59-A014 



Pin, Locating 

5 

11 59-DO 15 



Shim, Sealing Spacer 

1 

1 159-A016 



Ring, Retaining 

1 

Precision 

No. 

010-12137 

Packing "0" Ring 

3 

Precision 

No. 

161-12137 

Packing "0" Ring 

2 

Prec i s i on 

No. 

012-12137 

packing "0" Ring 

2 

Free i s i on 

No, 

168-12137 

Packing M 0 M Ring 

3 

P rec i s i on 

No. 

253-12137 

Packing “0" Ring 

1 

Waldes No. 

N5000-775-S 

Retaining Ring 

1 

Waldes No. 

6900 

P 1 i ers 

1 




Screw, Cap Socket Head, No. 10-32 x 
1 in Long Stainless Steel 

23 




Tube (0.0625 O.D. x 0.010 Wall), 
Stainless Steel 

As 

Requi red 

100-1-1 



Fitting (0.0625 Tube x 0.0625 Pipe), 
Swage 1 ok 

2 

100-1-2 



Fitting (0.0625 Tube x 0.125 Pipe), 
Swage lok 

2 


For ail other components of turbine, see turbine research package assembly, 
Drawing No. 699801 of Reference 2. 
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FIGURES 



Station Nomenclature 


0 

1 

2 

3 

31 

4 

5 


Scroll Inlet 
Nozzle Inlet 
Nozzle Exi t 

Rotor Inlet (First Blade Row Inlet) 
Intermediate Station (Second Blade Row Inlet) 
Rotor Exi t 
Diffuser Exit 


F 1 CURE 1 “ SCHEMATIC PICTURE OF RADIAL INFLOW TURBINE STAGE 



Static Efficiency, 7? Static Efficiency, V 

4sS • 7 ’ U &S 


70 



0.96 0.98 1.00 1.02 1.04 


Ratio of Rotor Inlet Diameter to Datum Diameter 



Ratio of Nozzle Vane Inlet Diameter to Datum Diameter 
FIGURE 2 - STATIC EFFICIENCY VARIATION WITH NOZZLE INLET DIAMETERS FOR 

three" 'datum' designs and variation of static efficiency, disk, and 

BENDING STRESSES W I TH ROTOR I N LET D 1 AMETERS FOR DATUM 1 DES I GN 


Stress x 1 0~^ (*kN/m 2 ) 




Static Efficiency, % _ Static Efficiency 
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FIGURE 3 - VARIATION OF STAT I C EFF I C I ENCY WITH OUTER DIAMETER 
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FIGURE 7 " DESIGN VELOCITY DIAGRAMS FOR THE VANEP NOZZLE AND THE 

ROTOR INLET STATION 



























6 - 



FIGURE 8 - DESIGN VELOCITY D-IAGRAMS AT THE ROTOR DISCHARGE STATION 
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FI GURE 9 - THE BLADE LOADING OBJECTIVE FOR THE OUTER CASING 
STREAMTUBE OF THE REDESIGNED TURBINE (16 BLADES IN EACH 
ROW OF A TANDEM' ROW' CONFIGURATION) 
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FIGURE 10 - FINAL MERIDIONAL FLOW PATH OF THE BRU TURBINE ROTOR 
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FIGURE 13 - FIRST BLADE ROW SURFACE VELOCITY DISTRIBUTION 
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FIGURE 14 - SECOND BLADE ROW SURFACE MACH NUMBER DISTRIBUTION AND 
CORRESPONDING BLADE SHAPE AT A RADIUS OF 1.70 INCHES 
(NEAR THE OUTER CONTOUR) “ ' “ 
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FIGURE 15 n SECOND BLADE ROW SURFACE MACH NUMBER DISTRIBUTION AND 
CORRESPONDING BLADE SHAPE AT A RADIUS OF 1.125 INCHES 
— — (NEAR THE HUB) ' ■ ; ' / 
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FIGURE 18 - BLADE ANGLE VARIATION ALONG THE STREAMLINE NEAR THE OUTER CASING 


CASING STREAMTUBE CENTERS 
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FIGURE 19 - MERIDIONAL FLOW PATH OF THE NOZZLE 
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FIGURE 20 - SURFACE VELOCITY DISTRIBUTION OF THE NOZZLE VANE 












FIGURE 21 - NOZZLE VANE GEOMETRY SELECTED FOR THE FINAL DESIGN 
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FIGURE 22 - TURBI NE INLET SCROLL GEOMETRY 
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FIGURE 23 - TWO-DIMENSIONAL BOUNDARY-LAYER PARAMETERS AND SECONDARY FLOW 
" ' PARAMETER Cp OF THE FIRST BLADE ROW (CAS 1 NG-TthFAMTU BE) 
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FIGURE 28 - TWO-DIMENSIONAL BOUNDARY-LAYER PARAMETERS AND SECONDARY FLOW 
PARAMETER Cp OF THE SECOND BLADE ROW (HUB STREAMTUBE) ~~ 
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FIGURE 30 - TWO-DIMENSIONAL BOUNDARY-LAYER PARAMETERS AND SECONDARY FLOW 
PARAMETER & FROM HUB TO CAS 1 NG OF THE SECOND 
*" BLADE ROW (PRESSURE SURFACE) * 
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FIGURE 32 - TWO-DIMENSIONAL BOUNDARY-LAYER PARAMETERS AND SECONDARY FLOW 
PARAMETER Jg FROM CASING TO HUB ALONG THE SUCTION SURFACE 

OF THE NOZZLE VANE 


CM 



Incompressible Shape Factor (H 




CM 


0 0.2 0.4 0.6 0.8 1.0 

Normalized Camber Line Distance 

FIGURE 33 "TWO-DIMENSIONAL BOUNDARY-LAYER PARAMETERS AND SECONDARY FLOW 
' ' PARAMETER FROM CAS I NG TO HUB ALONG THE PRESSURE’ ' ' 

SURFACE OF THE NOZZLE VANE 
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FIGURE 36 - ROTOR TEMPERATURE DISTRIBUTION AT 100 PER CENT SPEED 
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FIGURE 40 - RADIAL STRESS DISTRIBUTION AT' DESIGN TEMPERATURE 
OPERATING AT 1 00 PER CENT SPEED 
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FIGURE 41 -TANGENTIAL STRESS DISTRIBUTION AT DESIGN TEMPERATURE 

OPERATING AT 100 PER CENT SPEED. 
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FIGURE 43 - RADIAL VARIATION OF AXIAL DEFLECTION AT DESIGN AND 
ROOM TEMPERATURES (REFERENCE ATg g Of 
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FIGURE kS - NORMAL THICKNESS DISTRIBUTION ALONG THE HUB AND 
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FIGURE 49 - CAMPBELL DIAGRAM FOR THE ROTOR 



















































Turbine End 
Bearing No. I 


Rear 

Bearing No. 2 



(3) The moment transmitted through the wafers was assumed negligible. 
Turbine Wheel In essence, a pinned point was assumed at each end of the flexible 

' coupling. 

Mass 3.8679 (1-7 544 kg) (4) Rear beari ng spring rate assumed constant, K = 2.5 x 10 5 Ib/in 

|p 7.7269 lb-in 2 (22.61 N/cm ) * 

kgcm 2 ) (5) Bearing 2 takes thrust load in direction shown. 

| 4.5647 lb-in 2 (13-36 (6) Bearing loads calculated as follows: (lbs) (lbs) Thrust (lbs) 

kgcm 2 ) 11.1 13.2 75 

(49.4 N) (58.7 N) (333.6 N) 


FIGURE 51 - TURBINE ROTOR SYSTEM CHARACTERISTICS 
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FIGURE 53 - RADIAL VARIATION OF FRONTAL AND BACK PRESSURES 

ON THE TURBI NE WHEEL 
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NOMENCLATURE 


Symbols 

A 

b 

et 

C 4 

C w 

O p 

D 

D w 







Description 


Units 


Area 

Passage width 
Biot number 

Skin friction coefficient 

Meridional velocity 

Specific heat at constant 
pressure 

Diameter 

Hydraulic diameter 

Diffusion parameter (Equa- 
tion 2) 

Frequency 

Gravitational constant 
Heat transfer coefficient 
Incompressible shape factor, 

O 

Energy conversion factor 
Axial length 
Length 
Mass flow 

Rotational speed 
Specific speed, Nfq 

ahV 

Pressure 
Prandtl number 
Volume flow rate 
Reynolds number, 

Radius 

Blade spacing 


sq ft (sq cm) 
in (cm) 


fps (m/sec) 

Btu/lbm deg R 
i n (cm) 
in (cm) 


cycles per 
sec 

Ibm ft/lbf 
sec^ 

Btu/hr sq ft 
deg F 


ft lbf/Btu 


i n (cm) 

i n (cm) 

Ibm/sec (kg/ 
sec) 


rpm 

ft 3/4 


/(min) 

(sec)^'^ (m3/^/ 
(mi n) (sec) V2) 


psia 


(kN/m 2 ) 


cu ft/sec 


ft (cm) 
i n (cm) 
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Symbols 

T 

U 

V 

w 

E 

o< 

€ 

AH 

/i 

? 

CO 

0 

n, 

Cp 

Subscri pts 

b 

Cvtb 

a 

x 

ft 

7n 

n 

© 

ov 

Ref 

*h 


Descri ption 


Uni ts 


Temperature 

Trailing edge thickness 

Wheel speed 

Absolute velocity 

Relative velocity 

Number of blades 

Absolute flow angle from axial 
di recti on 

Displacement thickness 

Isentropic head on total pres- 
sure ratio 

Dynamic viscosity 
Density 

Angular velocity 
Momentum thickness 
Eff i c i ency 

Secondary flow parameter 
Secondary flow parameter 


deg R (deg K) 
i n (cm) 
fps (m/sec) 
fps (m/sec) 
fps (m/sec) 


deg 

i n (cm) 

ft (m) 

Ibm/ft sec 
(kg/m sec) 

lbm/ft^ (g/ 
citH) 

rad/sec 

i n (cm) 


» '<* «*/ 'L 6 ^ 

,4CW^-Wkv - 
■■ ■ ^ 


Description 


Back 

Conditions corresponding to 
Mach number of 1.0 

Driving surface 

Frontal 

Hub 

Mean 

Nozz 1 e 

Total condi t ions 
Over-al 1 
Reference 
Rotor 
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Subscripts 


O 

t 

a. 

3 
31 

4 

5 


Descrl pt ion 

Shroud 

Static 

Trai 1 i ng surface 
Wal 1 

Station at scroll inlet 
Station at nozzle inlet 
Station at nozzle exit 
Station at rotor inlet 
Intermediate station 
Station at rotor exit 
Station at diffuser exit 
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APPENDIX I 

AERODYNAMIC DESIGN PROCEDURES 
Stator Aerodynamic Design Method 

Scroll Design Procedure 

The method used for the detailed design of the scroll is based 
on Reference 17 and may be used for arbitrary cross-sectional shapes. The 
method determines the cross-sectional area of a scroll, as a function of 
the azimuth angle, that is necessary to obtain uniform static pressure 
distribution around the circumference. 

Assumptions and Limitations 

The method basically solves the continuity equation and angu- 
lar momentum equation with the following assumptions: 
t l . Flow is incompressible. 

2. Free vortex flow. 

3. FJow is one-dimensional. 

4. No secondary flow effects. 

The cross section of the scroll is calculated in two stages. 

1. Calculation of area at a specified azimuth angle with the 
above assumptions. 

2. Correction of the above area to allow for frictional ef- 
fects. 

The method is limited to cases of low Mach numbers with uniform flows. 

Input and Output Description 

The following are the input requirements in addition to the 
specification of type of cross section and any dimensional restrictions 
in the radial direction: 

1, Azimuth angle at scroll inlet. 

2. Stator inlet flow width. 
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3. Stator inlet flow angle. 

4. Stator Inlet radius. 

5. Friction factor. 

The output values at a given azimuth angle consist of: 

1. Mean radius of the cross section. 

2. Corrected area of the cross section. 

3. Height of the cross section. 

4. Head loss. 

Nozzle Aerodynamic Design Procedure 

The method employed in the aerodynamic design of the rotor ap- 
plies here in essence. An important simplification in stator design is 
introduced by the use of two-dimensional blades, which may be specified 
directly in terms of camber line angle instead of indirectly by prescribed 
loading. The method reduces to an analysis of the flow field and boundary 
layers for specified passage geometry, which is systematically modified 
until satisfactory performance is indicated. A description of the rele- 
vant design tools, their capabilities and limitations and the required 
input and resulting output are presented in the discussion of rotor design 
method . 


Vaneless Space Design Procedure 

The flow analysis of the vaneless space is based on one-dimensional 
compressible flow equations in an axi symmetric radial duct as described in 
Reference 18. The angular momentum and continuity equations are solved with 
the aid of an empirical expression for the friction coefficient as a func- 
tion of Reynolds and Mach numbers. An iterative integration of these 
equations has been chosen as the method most applicable to digital computer 
computations and a corresponding routine (Program VLS) is available. The 
chief limitation of this method is the assumption of one-dimensional flow. 

Input data are the axial channel width as a function of radius, 
the viscosity of the fluid and its variation with temperature, and inlet 
conditions consisting of the stagnation state, inlet Mach number, and flow 
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angle. The output comprises the flow conditions and the static and total 
state at specified radial stations. 


Rotor Aerodynamic Design Method 


Design Procedure 


The detailed aerodynamic design of the turbine rotor involves 
the fol lowing steps: 

1. Generation of a blade surface from prescribed loading dis- 
tribution and imposed geometrical constraints. 

2. Quasi -three-dimensional analysis of the inviscid flow field 
In the turbine blade passages. 

3* Detailed three-dimensional flow analysis in the blade row 
inlet and discharge regions. 

4. Compressible turbulent boundary- layer analysis on the walls 
of the blade passage and evaluation of discharge losses. 

The initial blade loading specifications may be related to the rotor per- 
formance and systematically revised by repeating the above series of steps 
in order to produce an optimum design. Concurrently, the geometric con- 
straints imposed on the design of the blade surface must be evaluated or 
modified in terms of permissible blade root bending stresses and manufac- 
turing feas ibi 1 i ty . 

NREC has developed computing routines for this design process 
which include the following: 

1. Program DESIGN to generate detailed blade geometry from 
given blade loading, flow path, and geometrical specifica- 
tions. The theoretical basis for this method has been re- 
ported in Reference 19 and the computer program in Refer- 
ence 20. 

2. Program ANALS to analyze the three-dimensional inviscid 
flow in a rotor blade passage of specified geometry. This 
method is also described in detail in References 19 and 20. 
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3. Program BL2BL to analyze In detail the inviscid flow field 
of the blade end regions. This method has been reported in 
References 21 and 22. 

4. Program BIANCA to compute compressible turbulent boundary- 
layer characteristics. This method is described in Refer- 
ence 23. 

Basis, Capabilities, and Limitations of Design Methods 

The basis of NREC's current Impeller blade design method is the 
assumption that the flow can be analyzed In ax i symmetric streamtubes. This 
amounts to ignoring the effects of viscid and inviscid secondary flow and 
represents a limitation. 

The determination of the blade shape is based on the relation 
between blade surface velocities, streamline shape, and streamwise dis- 
tribution of angular momentum (or loading) as expressed by the momentum 
equation. The prescribed distribution of loading along any streamline 
defines for a known meridional flow field, the geometry of that streamline. 
When corrections are applied to that geometry to account for incidence and 
deviation near the blade end regions, the corresponding camber line is 
found. This process may be repeated for one or two additional streamlines 
in the flow passage. A geometrical constraint is imposed on the loading 
of additional streamlines by the requirement that in order to define a 
practical blade shape, the several camber lines must be aligned on speci- 
fied leading or trailing edges. The remainder of the blade surface is 
then determined by the specification of straight-line generatrices. Such ’ 
straight-line elements are uniquely determined by three camber lines. An 
alternate blade element specification, which is frequently used, is to 
prescribe along one camber line the element orientation relative to the 
axis and the meridional plane. 

The limitations of this method are evidently the restriction 
to straight-line elements and certain approximations in the use of the 
momentum equation and in the assumed streamwise variation of deviation be- 
tween flow angle and blade camber line angle inside the blade passage. 
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The flow field analysis in the blade passages determines the 
positions of the axisymmetric streamtubes in the blade passage by appli- 
cation of the momentum, continuity, and energy equations. This result 
is achieved by an iterative procedure in which the positions of the 
streamlines which define the streamtubes of specified mass flow incre- 
ments are initially estimated and subsequently revised by the flow field 
solution corresponding to their geometric properties of slope and curva- 
ture. When successive streamline estimates agree to within a specified 
tolerance, the solution has converged. In a second step this axisymmetric 
solution is extended in circumferential direction to the blade driving and 
trailing surfaces by using approximations to the momentum equation first 
applied in Reference 27. 

The limitations of this method lie chiefly in the latter ap- 
proximations which preclude the determination of the deviation between 
mean flow angle and blade camber line angle. In the major part of the 
blade passage this deviation may be negligible. Near the blade end re- 
gion the computed blade surface conditions are certainly inaccurate. A 
thorough analysis of the flow field in these areas is accomplished by the 
following method which is designed to calculate the compressible flow 
field within an infinitesimal axisymmetric streamtube of variable thick- 
ness, as calculated by the foregoing method. The approach here is to 
solve the two-dimensional stream function equations by a relaxation 
method; NREC's computer program (BL2BL) is based on the method developed 
by T. Katsanis of NASA LeRC (Ref 5)* A limitation of the method is the 
restriction to subsonic flows. 

The analysis of the compressible turbulent boundary layers along 
the passages of the blade surface is based on an explicit integral theory, 
which is characterized by the initial integration of the governing equa- 
tions across the boundary layer and the computation of integrated parame- 
ters from the resulting ordinary differential equations. The essence of 
the method is the solution of the momentum and kinetic energy equations 
integrated across the boundary layer. These basic relations are supple- 
mented by semi -emp i ri cal relations for the wall shear stress coefficient 
and the dissipation coefficient. 
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Input Requirements 

The input information required for the blade design (Program 
DESIGN) consists of hub and shroud and blade edge contours in the 
meridional plane, the normal blade thicknesses along the hub and shroud, 
the blade element orientation at leading and trailing edges, inlet flow 
conditions consisting of total pressure, temperature and angular momen- 
tum, an average polytropic efficiency, the total mass flow, the operat- 
ing speed, and the specification of incidence and the angular momentum 
at discharge. 

The flow analysis Input data (Program ANALS) consists of the 
inlet flow conditions and operating conditions (as above) and the specifi- 
cation of the flow path and blade shape. The description of the blade 
surface is in the form of a grid of points defined by the cylindrical 
point coordinates and local tangential blade thickness. Additional data 
are the incidence factor (defined like the slip factor), the location 
along the blade camber lines where the flow deviation ceases, and a poly- 
tropic efficiency. 

Input data for detailed flow analysis in the blade end region 
(Program BL2BL) are streamtube and blade section geometry, fluid proper- 
ties and flow conditions, the distribution of losses along the flow path, 
and the spatial domain and finite difference mesh geometry. Flow condi- 
tions are specified by the absolute total temperature, total pressure, 
entropy and swirl upstream of the blade, and the streamtube mass flow or 
upstream flow angle. 

Input data for the compressible boundary-layer analysis (Pro- 
gram BIANCA) consists of streamline geometry on the passage surface, fluid 
properties, and initial and boundary conditions. The geometric informa- 
tion consists of the definition of the free streamlines along which 
boundary-layer properties are to be computed, including radius and ori- 
entation of the normal with respect to the axis of rotation as a function 
of streamline distance, and the width of the streamtubes (measured on the 
surface) for which the streamlines are representative. Fluid property in- 
formation required includes the gas constant, specific heat and specific 



132 


heat ratio, and fluid viscosity at standard temperature. Boundary condN 
tions required are, in general, the free stream velocity and static pres- 
sure distribution along the streamlines and the rotational speed of the 
surface. 

Output 

1. Full specification of flow path and blade geometry. 

2. Free stream flow conditions on the blade surfaces. 

3. Boundary-layer characteristics along representative stream- 
lines of the blade passage and passage discharge losses. 

Diffuser Aerodynamic Design Method 

Design Procedure 

The detailed design of annular diffusers of optimized over-all 
dimensions reduces to the proper selection of cross-sectional area as a 
function of diffuser length. The basis for this choice is an adequate 
description of the internal flow field and of the boundary layers. NREC 
has formulated a procedure to compute diffuser effectiveness and the 
point of separation for two-dimensional and annular diffusers. The 
method was designed to allow for the effect of a distorted velocity pro- 
file at inlet, but not for curvature terms and tangential swirl. A com- 
parison with experimental results shows that the method can predict dif- 
fuser effectiveness with fair accuracy for diffusers with uniform flow 
at inlet, and also predict the development of a shear profile, if the 
inlet profile is in equilibrium. 

The method involves a streamtube analysis to calculate the 
development of the velocity field. The output from this calculation is 
used to obtain the characteristics of the boundary layer, and hence the 
point of flow separation. The following steps may be identified: 

1. Calculate the stagnation properties of a specified number 
of streamtubes at diffuser inlet. 
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2. Make a first estimate of the boundary-layer displacement 
thickness to give the effective diffuser wail position. 

3. Proceed downstream In steps, calculating the pressure and 
velocity profile at each axial station. 

4. Using the pressure distribution thus obtained, make new 
estimates of the boundary- layer displacement thickness and, 
hence, the displaced wal 1 coordinates. 

5» Steps 3 and 4 are repeated until the solution converges. 

A computer routine is available to automate this procedure and 
is described in References 25 and 26. 

Step 4 in the above procedure merits further comment here. A 
key parameter is the shape factor, H j which is the ratio of the displace- 
ment thickness to the momentum thickness 

H- % 


Using the velocity distribution obtained in Step 3, the momentum thick- 
ness is obtained by solving a form of the momentum integral equation for 
compressible flow. A simple empirical relation suggested by Dussourd 
(Ref 27) from work on diffusers is used to find the shape factor. These 
two relations are presented in Appendix I of Reference 25. The new dis- 
placement thickness is then obtained directly. 

The separation criterion is expressed in terms of the shape 
factor, (-{ . Separation occurs when {-f reaches a critical value, f4sg{5 
Doenhoff and Tetervin (Ref 28) measured HggpOn NACA-65 aerofoils at mod- 
erate incidence angles and found values in the range 1.8 to 2.2. 

To use the streamtube method, a value of Hggp* 1*9 is recom- 
mended, This value was selected to give the best correlation between the 
streamtube method and the experimental results of Reneau, Johnston, and 
Kline (Ref 29) for two-dimensional diffusers, and Sovran and (Clomp- (Ref 
30) for annular diffusers. 


Basis, Assumptions, and Limitations 
1. There is no radial pressure gradient or tangential swirl. 
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2. The flow may be represented by a number of streamtubes. 

The flow in each streamtube is uniform and isentropic, and 

no mixing occurs between adjacent streamtubes. The fluid 
obeys the perfect gas law. 

3 . The boundary layer is identical to that on a flat plate 
with the same pressure gradient that exists in the diffuser. 
The effect of the boundary layer is to produce an effective 
displacement of the diffuser wall. When the flow separates, 
no further pressure rise occurs. 

4. There is a linear variation of streamline slope from one 

wall to the other at any axial station. At the first sta- 

tion there is zero streamline slope. 

A limitation of the procedure is that the inviscid flow calcu- 
lation does not allow for curvature terms and tangential swirl. 

Input Requirements 

The input conditions at the inlet station are: 

1. Weight mean total pressure. 

2. Uniform stagnation temperature. 

3 . Mass flow rate. 

4. A velocity profile given by 

= TOD 

where IA.- an arbitrary velocity 

a tabulated function of-^- sga t nst ^ . 

5 . The boundary-layer displacement thickness. 

6 . The boundary-layer shape factor. 

Output 

Output data at intermediate and final stations are: 

1. The velocity profile. 

2. The boundary-layer displacement thickness. 

3. The boundary- layer shape factor. 
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APPENDIX II 

MECHANICAL DESIGN PROCEDURES 


Rotor Mechanical Design Methods 

Rotor Steady-State Heat Transfer Coefficients 

In order to predict the temperature distribution in the rotor, 
it is necessary to compute the surface temperature of the metal. These 
temperatures depend upon the gas and adiabatic wall temperatures as well 
as on the heat transfer coefficients. Both heat transfer coefficients 
and Biot numbers are calculated. For the disk surface, the heat trans- 
fer coefficient may be calculated by standard methods. It will in general 
depend upon the Reynolds and Prandtl numbers and is computed by equations 
derived empirically for high velocity flow. For the areas where blades 
exist, the surface is treated as a finned surface. The additional heat 
transfer is accounted for by correction of the standard heat transfer 
coefficients derived above. 


Limi tat ions 

The major limitation of this process is in the selection of the 
formula for calculation of the heat transfer coefficient. The experimental 
conditions leading to the equation should be as similar to flow in the im- 
pel ler as possible. 


Input Requirements 

1 . Impel ler geometry. 

2. Gas properties. 

3. Gas temperatures, velocity, and pressure distribution. 

Output 

1. Heat transfer coefficient distribution. 
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2. Biot number distribution. 

3« Adiabatic wal 1 temperature distribution. 

Calculation of Temperature Distribution in Rotor 

The calculation of the rotor temperature distribution is accom- 
plished when the Biot number and adiabatic wall temperature distribution 
are known. It is desired to determine the metal surface temperature dis- 
tribution. This is accomplished by NREC Program STRESS (Refs 3 and 4) 
which uses an "energy-integral" temperature analysis procedure. The im- 
peller is divided into a series of concentric cylinders for the analysis. 
A typical annular cylindrical element will have heat flow in or out of 
Its four surfaces. All elements are then solved simultaneously with the 
boundary conditions included. 

Limitations 

The axial temperature distribution is an approximation and its 
accuracy depends upon the selection of appropriate constants. The method 
also depends upon the Initial calculation of Biot numbers and temperature 
distributions. 


input Requirements 

1. Biot number distribution. 

2. Adiabatic wall temperature distribution. 

3« Rotor geometry. 

4. Rotor material properties. 

Output 

Rotor surface temperature distribution. 

Thermal and Mechanical Disk Stress 
Impeller disk stresses are calculated using NREC Program STRESS 
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(Ref 4). Both thermal and centrifugal loadings are considered. The basic 
method employed in the analysis is that of Schilhansl. To use the method, 
the rotor geometry, rotor face temperature distribution, and material 
properties must be known. The computer program uses a series of concen- 
tric cylinders to determine centrifugal and thermal distributions and 
loadings. The disk back and front face stress as well as radial and axial 
disk deformations are computed. 

For burst considerations the average tangential stress at the 
worst location will be computed by an area integration method. 

Limi tations 

The approximation itself and the assumption of linear cylindri- 
cal deformations is the major limitation of the method. The procedure is 
thus not good for the axial portions of long rotors. The results will 
show a low stress due to centrifugal loading at the impeller eye. 

Input Requirements 

1. Impeller geometry and operating conditions. 

2. Impeller face temperature distribution. 

3- Material properties. 

Output 

1. Radial and tangential disk stresses for the impeller back 
and front face. 

2. Radial d i sk deformat i on at the axial plane of the turbine 
Inlet. 

3. Axial disk deformation at the axial plane of the turbine 
inlet. 

4. Average tangential stress at 140 per cent speed. 

Blade Root Stresses 

To predict the blade root stress due to centrifugal loading, 
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NREC Program BLDROT (Ref 14) is used. The method employed by the program 
Is a rectangular beam approximation. That is, the continuous blade is 
approximated by a series of beams of rectangular cross section. The beams 
are pinned at the shroud line to adjacent beams in a direction normal to 
the blade. Motions other than normal are thus unrestrained and lead to 
conservative results. 


Limitations 

The beam approximation method itself is a limitation. It is, 
however, not a severe limitation as it has been found that the magnitude 
of the stress computed is as would be expected. The location of the 
stress, however, particularly near the boundaries is not accurate. 

Another limitation is that the stresses calculated by the pro- 
gram are at present limited to bending and tensile. This is generally 
sufficient, particularly where nonradial blade elements are considered 
as these stresses will predominate. In radial blades, if radial beam 
elements are used, a conservative estimate of stress will be obtained 
since the radial direction would be unrestrained. 

Input Requirements 

1. Geometric definition of blades. 

2. Operating speeds. 

3- Material properties. 

Output 

1. Blade root tensile stress distribution. 

2, Blade root bending stress distribution. 

Blade Natural Frequency 

Impeller blade natural frequencies, in general, are much higher 
than operating speeds. Where this is the case, relatively crude methods 
may be used to estimate natural frequencies. The methods used employ 
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plate theory. A section of the blade may be thought of as a plate, 
clamped on two sides which would represent its attachment to the hub and 
the adjacent blade section. The eye of the impeller is generally the area 
of lowest frequency. This is treated separately by a method due to Voysey 
(Ref 9). This method is also a plate oriented method but considers taper 
in a more rigorous fashion and has been checked experimentally. It includes 
the experimental factor derived to bring the analytical results into agree- 
ment with the empirical. 


Limitations 

The major limitation of the method discussed is that segments 
of the blade are treated separately. In this manner, interaction effects 
with the remainder of the blade are either ignored or approximated. 

If a blade frequency is in the range of an operating or funda- 
mental frequency, it will be necessary to perform a more sophisticated 
analysis using finite element techniques. In this manner a continuous 
blade analysis may be performed. 

Input Requirements 

1. Blade geometry. 

2. Material property. 

3* Operating speed. 


Output 

Impeller blade first natural frequency. 

Blade Vibratory Stress Levels 

The vibratory blade stresses, or those caused by pulsating 
aerodynamic loads are, in general, quite small when compared with the 
centrifugal loading. They do become important, however, when operation 
occurs near a natural blade resonance. To estimate these stresses, it 
is necessary to know actual blade aerodynamic loading. This may be 
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doubled since the worst variation will be from zero to twice the steady- 
state loading. To determine an amplification factor the method of Trumpler 
and Owens (Ref 15) may be used. They give an amplification factor to be 
applied to the stress result based on a comparison of the frequency of the 
disturbance with that of the fundamental mode. 

In general, the target value for maximum vibratory stress is 
10,000 psi or less. 


Limitations 

The limitation of the method is that the amplification factor 
is a statistically derived number based on experimental results. A lot 
depends upon the manufacturing quality; that is, the existence of stress 
risers and the like. In addition, the fundamental frequencies themselves 
are not that well known. These limitations are not significant where 
resonance and operation are widely separated. 

Input Requirements 

1. Gas blade loading magnitude. 

2. Natural frequencies of interest of the blade. 

3. Disturbance frequencies. 

Output 

Estimate of vibratory stresses at the blade root. 

Flexural Vibrations (Critical Speeds) 

Flexural vibrations in the nature of critical speed and rela- 
tive mode shapes will be computed for the turbine rotor system. The sys- 
tem to be considered consists of the rotor and shaft in two bearings. 

The method employed is a modified '’Holzer 11 technique (Ref 16) using the 
NREC computer Program HOLZER. 

The program uses a stiffness approach using f ini te elements . 

The finite elements are lump masses joined to massless bars. The rotor 
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assembly Is thus approximated by a system of masses and bars. The bars 
will support both bending and shear. The lump masses may contain both 
mass and gyroscopic inertia. The two bearing stations may contain pre- 
specified spring gradients, a different gradient at each bearing. The 
cross section of the massless bars is arbitrary, since only moment of 
inertia need be known. 


LImi tations 

One limitation of the method is that it is basically an approxi- 
mation. Lump masses are used to approximate a distributed mass system. 
This is not considered a severe limitation as the approximation is quite 
good. Another limitation is the fact that damping at the bearings is not 
considered. Depending upon the bearing system, this may or may not be a 
severe limitation. In the bearing system considered for this application, 
damping is low such that the limitation is not severe. 

Input Requirements 

1. Approximate system dimensions. 

2. Masses. 

3. Shear stiffness. 

4. Bending stiffness. 

5. Impeller gyroscopic inertia. 

Output 

1. Natural frequencies for a practical range of impeller end 
spring constants for the final shaft impeller geometry. 

Mode shapes for the lower frequencies having the greatest 
effect on the operating speeds. 


2 . 
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Shaft Attachment Methods 


Design Procedure 

The shaft attachment method considered at present is a press 
fit. The impeller will have a stub shaft and the shaft will have a 
cavity to receive it. Those parameters needing calculation are the 
torque transmission capability under centrifugal load and the impeller 
and shaft stress under no external load. The procedure to be followed 
is detailed in Reference 31. 


Limi tations 

The major limitation in calculations such as this is that of 
accurately estimating the friction coefficient for torque transmission. 
This can be overcome by using a reasonably conservative estimate for 
friction coefficients in addition to insuring that the joint is not over- 
stressed. A reasonable target maximum figure is 3/4 yield strength of 
the material as the maximum joint stress. 

Input Requirements 

1. Torque transmission requirements. 

2. Shaft size. 

3. Material stress limits. 

Output 

1. The design and drawings of the impeller shaft assembly. 

2. Results of calculations for torque transmission capability 
of the joint. 

3'. Results of calculations for static joint stress conditions. 

Bearing Selection and Analysis 

Prior to ball bearing selection and analysis, a calculation of 
the anticipated thrust load variation will be made. The thrust loads are 
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computed by integrating the expected static pressure distribution along 
the front face and back face of the impeller to determine the resultant 
thrust force* On the back face, any anticipated pumping is considered. 
Having the expected thrust and radial loads, the bearings may be selected 
on the basis of Bjq life expectancies. The method of calculating the 
bearing life estimate is based on the method established by the Anti- 
Friction Bearing Manufacturers Association (AFBMA) . After these calcu- 
lations have been made, there is to be a comparison made with the 
existing bearing set. If major differences in Bjq life are found and 
it is possible to do so, the bearings will be reselected. 

Limitations 

The only major limitation is that in computing Bjq life a 
small percentage of the bearings fail, perhaps well in advance of the 
estimated life. In addition, care must be taken in assembly with regard 
to such factors as alignment, preload, and lubrication system to ensure 
that the calculations are valid. 

Input Requirements 

1, Results of critical speed calculations. 

2, Design data for existing system. 

3. Static pressure distributions at maximum and minimum ex- 
pected thrust loads. 

4 . Dynamic and static radial loads. 

Output 

1. A specification for the bearings. 

2. Rotor thrust load predictions. 

3. B, 0 1 ife predictions for the bearings. 
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Turbine Scroll, Nozzle, and Diffuser Design Methods 

Design Procedure 

The mechanical design of the stationary parts of the stage in- 
volved the calculation of the following items as a minimum requirement 
for a satisfactory end product: 

1. Surface temperatures at certain critical locations from 
the inlet to the outlet flanges. 

2, Stress and/or deflection at corresponding critical loca- 
tions from the inlet to the outlet flanges, as a function 
of local temperature. 

All temper- tures will be based upon steady-state conditions. 

The procedures use^ ur stationary parts such as the inlet scroll, the 
nozzle row, and the exhaust diffuser are relatively standard practice. 

The detailed methods of analysis are largely based on the standard text- 
books of References 8 and 32. In the present design, the stationary 
parts have been designed for the turbine research package rather than the 
BRU operating at high temperatures and high pressure differentials. Hence, 
relatively little effort is needed to be devoted to the mechanical analyses 
of the stationary parts. Rather, the emphasis was placed on designing a 
research vehicle which provided ease of assembly and modification of de- 
sign standard. 
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APPENDIX HI 

REPRODUCTIONS OF THE MANUFACTURING DRAWINGS 

The following pages contain reproductions of the drawings 
produced for the new turbine stage. These drawings have been reduced 
by arbitrarily selected scaling factors for inclusion in this report. 
The following Is a list of drawing numbers, drawing titles, and the 
page on which the drawing is reproduced. 


Drawing No. 

Title 

Page No. 

1 159-R001 

Research Package Layout 

146 

1 159-R002 

Research Package Assembly, Turbine 

147 

1159-R003 

Research Package Outline, Turbine 

148 

1159-D004 

Wheel Turbine 

149 

1159-E005 

Wheel Assembly, Turbine 

151 

1159-RQ06 

Scroll Fabrication Assembly 

152 

1 1 59-R007 

Scroll Machining Assembly 

153 

1159-R008 

Nozzle Maching Assembly 

154 

1 159-0009 

Shroud Machining 

155 

11 59-DO 10 

Diffuser Machining 

156 

11 59-DO 11 

Exhaust Adapter 

157 

1 159-0012 

Rotor Back Shield 

15$ 

11 59-BO 13 

Pin Cap 

159 

1159-A014 

Pin, Locating 

160 

1159-DO 1 5 

Shim, Seal Ing Spacer 

161 

1 1 59-A01 6 

Retaining Ring 

162 
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